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Abstract 

This  research  deals  with  the  theoretical  prediction  of  the  overall  or  effective  behavior  of 
nonlinear  composite  materials  undergoing  large  deformations.  In  particular,  appUcations  are 
envisaged  to  the  high-temperature  creeping  behavior  of  metal/metal,  metal/ceramic 
composites,  and  to  porous  materials.  Both  isotropic  and  anisotropic  configurations  are 
considered  including  the  technologically  important  case  of  fibre-reinforced  composites.  The 
approach  is  based  on  new  variational  principles  developed  recently  by  the  author  (under 
AFOSR  sponsorship),  which  allow  the  estimation  of  the  overall  behavior  of  a  given 
nonlinear  composite  in  terms  of  the  effective  properties  of  a  suitably  optimized  linear 
comparison  composite  (with  the  same  microstructure).  The  key  advantage  of  the  method  is 
that  it  allows  direct  application  of  the  extensive  literature  on  linear  composite  materials,  in 
the  form  of  estimates  of  various  types  and  rigorous  bounds,  to  obtain  corresponding  results 
for  nonlinear  composites.  Additionally,  the  procedure  is  remarkably  simple  to  implement, 
and  the  final  results  are  usually  expressed  in  terms  of  finite-optimization  problems,  which 
can  be  readily  solved  with  little  compi  ational  effort.  Recent  progress  includes  the 
application  of  the  method  to  the  computation  of  the  effective  behavior  of  anisotropic 
systems,  including  laminated  and  fiber-reinforced  nonlinear  composites,  and  to  the 
determination  of  extremal  yield  surfaces  for  rigid/plastic  systems. 


Research  goals 

The  central  goal  of  this  project  is  to  estimate  the  effective  constitutive  properties  of 
nonlinear  composite  materials  undergoing  large  deformations.  Two  types  of  large 
deformations  are  of  particular  interest:  large  elastic  deformations,  corresponding  to 
materials  such  as  polymeric  composites,  rubber  foams  and  solid  rocket  fuel  composites; 
and  large  viscous  deformations,  corresponding  to  the  high-temperature  creeping,  or  to  the 
dynamic  plastic  deformation  of  metals.  The  main  current  goal — and  the  goal  for  the  next 
funding  period — will  be  to  estimate  the  effective  behavior  of  nonlinear  anisotropic 
composites,  such  as  fiber-re'  aforced/metal-matrix  composites.  Preliminary  results  (mostly 
for  laminated  materials)  indicate  that  there  may  be  significant  coupling  between  different 
loading  modes. 
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BACKGROUND 


There  is  a  fairly  recent,  but  growing  literature  in  the  field  of  composite  materials  that  is 
concerned  with  the  prediction  of  the  effective  behavior  of  nonlinear  composites.  To  date 
most  of  the  work  on  the  subject  has  made  use  of  ad  hoc  models  for  certain  classes  of 
composites,  as  well  as  of  heavy  numerical  computations  for  composites  with  periodic 
microstructures.  For  example,  Duva  and  Hutchinson  [7]  determined  dilute  estimates  for 
porous  creeping  composites  by  making  use  of  the  results  of  Budiansky,  Hutchinson  and 
Slutsky  [3]  for  a  single  void  in  an  infinite  matrix  of  power-law  viscous  material.  Duva  [6] 
proposed  a  differential  self-consistent  estimate  to  model  the  reinforcing  effect  of  rigid 
particles  in  a  power-law  creeping  matrix.  More  recently,  Bao,  Hutchinson  and  McMeeking 
[2]  proposed  generalized  self-consistent  estimates  for  the  yield  stress  of  two-phase, 
rigid/perfectly  plastic  composites.  On  the  other  hand,  Christman,  Needleman  and  Suresh 
[4]  and  Bao,  Hutchinson  and  McMeeking  [1]  carried  out  computations  on  the  unit  cell  of  a 
periodic  microstructure  to  model  the  effect  of  particle  size  and  shape  on  the  effective 
behavior  of  rigidly  reinforced  power-law,  metal-matrix  composites. 

As  far  as  rigorous  methods  for  estimating  the  effective  behavior  of  composites  with 
random  microstructures  (real  composite  materials  are  hardly  periodic),  the  only  method 
available  until  recently  was  the  Talbot  and  Willis  [15]  extension  of  the  Hashin-Shtrikman 
variational  principles  to  composites  with  nonlinear  behavior.  Ponte  Castaneda  and  Willis 
[14]  applied  this  procedure  to  nonlinearly  viscous  materials,  and  gave  the  first  rigorous 
bounds  and  self-consistent  estimates  for  the  effective  properties  of  composites  in  power- 
law  creep.  Later,  Ponte  Castaneda  [8]  developed  corresponding  bounds  in  the  context  of 
finite  elasticity.  Very  recently,  under  the  sponsorship  of  AFOSR,  Ponte  Castaneda  [9] 
developed  a  new  general  method  for  estimating  and  bounding  the  effective  behavior  of 
nonlinear  composites.  The  method  expresses  the  effective  behavior  of  nonlinear  composites 
in  trms  of  an  optimization  problem  involving  the  effective  behavior  of  linear  composites 
materials  with  similar  microstructures.  This  is  very  useful  because  of  the  large  body  of 
results  in  the  form  of  bounds  and  estimates  of  different  types  that  are  currently  available  for 
lineart  composites.  This  method  and  some  of  its  applications  to  specific  types  of 
composites  are  discussed  in  the  following  section. 


ACCOMPLISHMENTS 


Introduction 

In  this  section,  we  give  a  brief  account  of  the  theoretical  foundations  of  the  new  variational 
method  developed  to  predict  the  effective  behavior  of  nonlinear  composites.  Explicit  results 
for  composites  with  overall  isotropy  are  given  in  references  [9,  10],  and  corresponding 
results  for  fiber-reinforced  and  laminated  microstructures  in  [1 1]  and  [5, 12],  respectively. 
In  particular,  we  note  that  our  results  for  simple  two-phase,  power-law  creep  composites 
compare  very  favorably  with  the  results  of  other  authors  obtained  by  means  of 
sophisticated  large-scale  computations  for  periodic  composites.  Our  results  also  have  the 
advantage  that  they  can  often  be  given  simple  analytical  forms.  For  example,  we  have  been 
able  to  obtain  an  improved  characterization  of  porous  isotropic  and  transversely  isotropic 
materials;  in  fact,  because  of  the  bounding  properties  of  our  results,  we  are  able  to  show 
that  the  commonly  used  Gurson  model  is  not  appropriate  outside  the  range  of  dilute 
porosity  and  high-triaxiality  for  which  it  was  originally  developed  [13].  Similarly,  we  have 
been  able  to  obtain  yield  surfaces  for  isotropic  and  transversely  isotropic  composites  with 
rigid/perfectly  plastic  phases.  For  the  benefit  of  the  reader,  some  of  the  above  mentioned 
references  [5, 10, 1 1, 12]  are  included  in  the  Appendix. 

The  method 
Effective  properties 

Consider  a  two-phase  heterogeneous  material  occupying  a  region  in  space  of  unit  volume 
12.  We  characterize  the  local  stress  potential  ^(a,x)  of  the  material  in  terms  of  the 
homogeneous  phase  potentials  (t>^'\o)  via  the  relation 

(1) 

r-l 

where  the  (r  =  1, 2)  are  the  characteristic  functions  of  the  two  phases.  We  will  assume 
that  the  two  phases  are  nonlinear,  incompressible  and  isotropic,  so  that  the  potentials 
depend  only  on  the  effective  stress  O’,  =  -v/yS-S,  where  S  is  the  deviator  of  a. 
Thus,  we  write 


<2) 


where  the  scalar-valued  functions  are  assumed  to  satisfy  some  rather  weak  hypothesis. 

For  exanq>le,  may  be  chosen  to  have  the  power-law  form 


fW  = 


1  . 
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(3) 


with  n,  o;  and  denoting  material  constants  with  different  values  in  each  phase. 
The  local  constitutive  relation  for  the  creeping  material  is  given  by 


e  =  |^(CT,x),  (4) 

da 

where  c  is  the  rate-of-deformation  (strain-rate)  tensor. 

Letting  X  and  E  denote  the  average  values  of  the  stress  and  strain-rate  fields, 
respectively,  the  effective,  or  overall,  behavior  of  the  heterogeneous  creeping  material  is 
given  by  the  relation 


E 


cfL' 


(5) 


where  the  effective  energy  of  the  composite  d>  follows  firom  the  principle  of  minimum 
complementary  energy,  which  can  be  stated  in  the  form 


<I)(Z)=  nun  f0(CT,x)dx,  (6) 

aeS(£)Ja 

where  5(Z)  =  {a\Gjj  j  =  0  in  Q,  and  G^jiij  =  on  is  the  set  of  statically  admissible 
stresses. 

In  this  work,  we  are  interested  in  predicting  the  effective  behavior  of  two-phase 
heterogeneous  materials,  as  described  above.  We  further  restrict  our  consideration  to 
composites  with  phases  in  given  volume  fractions 

^  f  (7) 


New  variational  principles 

The  new  variational  principles  [9, 1 1]  are  based  on  a  representation  of  the  potential  for 
the  nonlinear  isotropic  material  in  terms  of  the  potentials  of  a  family  of  linear  comparison 
materials.  This  representation  is  obtained  with  the  help  of  the  Legendre  transformation. 
Thus,  for  a  homogeneous  nonlinear  material  with  "superquadratic"  growth  in  its  potential 
(j) ,  and  certain  additional  convexity  hypothesis,  we  have  that 


(8) 


0(a,x)  =  max{^,(CT,x)  -  V(/f,x)}, 

li>0 

where 

V(;i,x)  =  max{0^(CT,x)-0(a,x)}  (9) 

<3 

and  where  is  the  potential  of  a  linear  comparison  material  with  shear  modulus  /t .  Note 
that  in  the  above  optimizations,  the  position  vector  x  is  fixed. 

The  new  variational  principle  is  obtained  by  averaging  the  relation  (8)  and  the 
complementary  energy  principle  (6) ,  and  can  be  expressed  in  the  form 


where 


<I>(X)  =  max 

p(x)20 


<i>,(2:)-Jy(/t.x)dx|, 


(10) 


<I>,(Z)=  nun  f^„(CT,x)dx  (11) 

aeS(X)J 

a 

is  the  effective  potential  of  a  linear  heterogeneous  comparison  material  with  local  potential 
and  arbitrarily  variable  shear  moduli  /i(x) . 

The  variational  principle  described  by  (10)  and  (1 1)  can  be  given  the  interpretation  of 
first  solving  a  linear  problem  for  an  arbitrarily  heterogeneous  linear  material,  and  then 
optimizing  with  respect  to  the  variations  in  moduli  to  account  for  the  nonlinearity  in  the 
actual  material.  This  suggests  that  if  the  fields  happen  to  be  constant  over  the  nonlinear 
phases  in  the  actual  composite,  then  the  variable  moduli  /i(x)  can  be  replaced  by  constant 
moduli  over  each  phase,  i.e.  by  and  in  phase  1  and  phase  2,  respectively.  More 
generally,  however,  we  have  the  following  lower  bound  for  d>,  namely, 

d)_(Z)=  max  {<D  (12) 

^(1)  p(2)^0 

where  and  correspond  to  relation  (9)  evaluated  for  each  of  the  homogeneous 

phase  potentials  and  <p^^\  and  <I>„  is  the  effective  energy  of  an  isotropic  two-phase 

linear  composite  with  phase  moduli  and  in  volume  fractions  c”’and 
respectively.  This  latter  expression  has  the  advantage  over  the  previous  more  general 
relation  that  it  involves  only  a  finite-dimensional  optimization  problem,  and  can  be  further 
used  in  conjunction  with  bounds  and  estimates  for  the  linear  composite  to  induce 
corresponding  results  for  the  nonlinear  composite.  We  illustrate  these  possibilities  in  the 
next  section  for  a  general  two-phase  incompressible  composite. 


The  estimate  (12)  of  the  previous  section  for  the  effective  energy  of  the  nonlinear 
incompressible  composite  can  be  used  in  conjunction  with  the  well-known  Hashin- 
Shtrikman  bounds  for  the  effective  properties  of  the  linear  comparison  composite  to  yield 
estimates  for  the  effective  potential  of  the  nonlinear  composite.  It  is  shown  in  reference  [11] 
that  the  result  of  this  calculation  takes  on  the  simple  form 

=  min{c<V‘'(5“')  +  c<"y"'(5<"')},  (13) 

where  and  characterize  the  phase  behaviors,  and  5“’  and  are  simple 
expressions  of  o)  involving  the  phase  volume  fractions  and  and  the  overall 
effective  stress  applied  to  the  composite  Z,.  For  instance,  if  and  the  Hashin- 

Shtrikman  lower  bound  is  used  for  d)^,  the  corresponding  expressions  for  5^*^  and 
become 

=.j(l-c^^'<u)'  +|c'V  I,  and  =  (H-c"W)Z,  (14) 

Depending  on  the  relative  strengths  of  the  two  nonlinear  potentials  and 
the  result  of  putting  (14)  into  (13)  would  either  lead  to  a  rigorous  lower  bound  for  the 
nonlinear  potential  d>,  or  alternatively,  to  an  upper  "estimate"  for  d>.  Other  examples  are 
detailed  in  the  Appendices. 

RESEARCH  PLANS  FOR  THE  FUTURE 

In  the  remaining  nine  months  of  support  of  this  project,  we  plan  to  complete  the 
development  of  effective  constitutive  relations  for  anisotropic  nonlinear  composites, 
including  fiber-reinforced  solids,  for  a  fairly  broad  range  of  constitutive  behaviors  for  the 
matrix  and  inclusion  phases.  We  also  plan  to  get  started  in  the  study  of  damage 
localizations  in  composites,  which  may  arise  because  of  the  finite  kinematics,  or 
alternatively  through  thermomechanical,  or  other  micro-mechanisms.  These  are  important 
in  the  development  of  realistic  models  for  the  behavior  of  composite  materials,  as  discussed 
in  the  recent  "Panel  Discussion  on  Homogeneization  and  Damage  Localization  in 
Heterogeneous  Media"  in  this  year's  Contractors  Meeting  on  Mechanics  of  Materials 
(Dayton,  October  1991).  Thus,  we  hope  to  set  the  stage  for  a  new  project  under  AFOSR 
sponsorship.  Other  potential  avenues  of  research  include  the  consideration  of  isotropic 


composites  with  nonlinear  anisotropic  phases  (i.e.  polycrystalline  materials),  and  studies 
of  texture  evolution. 

Concerning  the  development  of  human  resources,  we  expect  that  the  student  supported 
over  the  past  two  years  under  this  project,  Mr.  Gal  de  Botton,  will  complete  the 
requirements  for  his  Ph.D.  degreee  in  the  Department  of  Mechanical  Engineering  of  the 
University  of  Pennsylvania  sometime  towards  the  second  half  of  the  next  year. 
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THEIR  APPLICATION  TO  COMPOSITES  MATERIALS 
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Abstract 

A  new  variational  method  for  bounding  the  effective  properties  of  nonlinear  composites  with 
isotropic  phases,  proposed  by  PONTE  CASTANEDA  {J.  Mech.  Phys.  Solids  39,  45,  1991),  is 
given  here  full  variational  principle  status.  Two  dual  versions  of  the  new  variational  principle  are 
presented  and  their  equivalence  to  each  other,  and  to  the  classical  variational  principles,  is 
demonstrated.  The  variational  principles  are  used  to  determine  bounds  and  estimates  for  the 
effective  energy  functions  of  nonlinear  composites  with  prescribed  volume  fractions  in  the  context 
of  plasticity.  The  classical  bounds  of  Voigt  and  Reuss  for  completely  anisotropic  composites  are 
recovered  from  the  new  variational  principle.s  and  are  given  alternative,  simpler  forms.  Also,  use  of 
a  new  identity  allows  the  determination  of  novel,  simpler  forms  for  nonlinear  Hashin-Shtrikman 
bounds  and  estimates  for  isotropic,  and  transversely  isotropic,  fiber-reinforced  composites. 
Additionally,  third-order  bounds  of  the  Beran  type  are  determined  for  the  first  time  for  nonlinear 
composites.  The  question  of  optimality  of  these  bounds  is  discussed  briefly. 
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1.  Introduction 

A  new  variational  procedure  for  bounding  the  effective  properties  of  nonlinear  composites  with 
isotropic  phases  has  been  proposed  recently  by  PONTE  CASTANEDA  (1991a).  The  procedure 
allows  the  estimation  of  the  effective  properties  of  nonlinear  composites  in  terms  of  the  effective 
properties  of  an  appropriately  optimized,  linear  comparison  composite  with  the  same  underlying 
microstructures  as  the  nonlinear  composite.  The  general  nonlinear  estimate  can  be  expressed  in  the 
form  of  a  lower  bound  for  the  effective  complementary-energy  function  of  the  composite,  in  such  a 
way  that  if  a  lower  bound  is  available  for  the  linear  comparison  composite,  a  lower  bound  can  be 
deduced  for  the  nonlinear  composite.  On  the  other  hand,  an  upper  bound  for  the  effective 
complementary-energy  function  of  the  comparison  composite  does  not  yield  in  general  an  upper 
bound  for  the  nonlinear  composite.  In  this  sense,  the  new  procedure  is  similar  to  a  nonlinear 
generalization  of  the  Hashin-Shtrikman  procedure,  developed  by  TALBOT  and  WILLIS  (1985) 
following  WILLIS  (1983),  which  makes  use  of  a  linear  homogeneous  comparison  material  to 
obtain  bounds  and  estimates  of  the  Hashin-Shtrikman  type  for  nonlinear  composites.  The  new 
variational  procedure,  however,  allows  the  use  of  any  linear  bound,  or  estimate  for  the  comparison 
composite  to  yield  a  corresponding  nonlinear  bound,  or  estimate  for  the  nonlinear  composite. 

In  this  paper,  we  strengthen  the  variational  procedure  proposed  by  PONTE  CasTAI^DA 
(1991a)  to  a  full  variational  principle,  which  under  appropriate  hypothesis  can  be  shown  to  be 
completely  equivalent  to  the  classical  minimum  principles  of  potential  and  complementary  energy. 
This  goal  is  accomplished  essentially  by  allowing  the  comparison  composite  to  have  continuously 
varying  microstructure,  instead  of  being  constrained  to  have  the  same  (piecewise  constant) 
microstructure  as  the  nonlinear  composite.  Naturally,  the  v  itional  procedure  of  PONTE 
Castaneda  (1991a)  results  as  an  approximation,  but  more  generally,  the  new  variational 
principles  allows  the  computation  of  the  exact  properties  of  certain  classes  of  special 
microstructures.  For  example,  DE  BOTTON  and  PONTE  CASTANEDA  (1991)  have  made  use  of  the 
new  variational  principles  to  determine  simple  expressions  for  the  effective  properties  of  nonlinear 
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laminates.  Here,  we  will  make  use  of  the  exact  version  of  the  new  variational  principles  to 
demonstrate  the  optimality  of  some  of  the  nonlinear  bounds  developed  by  means  of  the 
approximate  version  of  the  variational  principle. 

The  paper  is  organized  as  follows.  First,  the  definition  of  the  effective  properties  of  nonlinear 
composites  is  reviewed.  Then,  two  dual  versions  of  the  new  variational  principle  are  presented, 
with  rigorous  proofs  of  their  respective  equivalence  to  each  other,  and  to  the  classical  variational 
principles.  In  the  following  section,  we  apply  the  approximate  version  of  the  new  variational 
principle,  together  with  a  new  identity  introduced  in  DE  BOTTON  and  PONTE  CASTAfiEDA  (1991), 
to  develop  remarkably  simple  forms  for  the  Voigt  and  Hashin-Shtrikman  lower  bounds  for  the 
effective  complementary-energy  function  of  nonlinear  composites  with  phases  in  prescribed 
volume  fractions.  Results  are  given  for  multiple-phase,  compressible  composites,  but  a  particularly 
simple  form  for  the  effective  energy  of  the  nonlinear  composite  is  obtained  for  the  case  of  two- 
phase  incompressible  composites.  These  results  are  expressed  in  terms  of  a  one-dimensional 
optimizations  problem  [see  equation  (4.34)],  and  generalize  the  findings  of  PONTE  CastaI^DA 
(1991b)  for  two-phase,  incompressible  composites  with  a  linear  and  a  nonlinear  phase.  For  the 
Hashin-Shtrikman  bounds,  results  are  given  not  only  for  isotropic  nonlinear  composites,  but  also 
for  transversely  isotropic,  fiber-reinforced  materials  [see  equation  (4.46)].  Additionally,  a  third- 
order  bound  of  the  Beran  type  is  given  for  isotropic  composites  [see  equation  (4.41)].  To  the 
knowledge  of  the  author,  this  is  the  first  higher-order  bound  ever  produced  in  the  context  of 
nonlinear  composites,  and  illustrates  the  versatility  of  the  new  variational  procedure.  Additionally, 
in  this  section,  we  give  a  derivation  of  the  Reuss  upper  bound  for  the  effective  complementary- 
energy  function  of  the  nonlinear  composite  by  means  of  the  exact  version  of  the  new  variational 
principle.  This  serves  to  demonstrate  that  the  exact  version  of  the  variational  principle  is  able  to 
deliver  upper  bounds  (which  cannot  be  obtained  from  the  approximate  version  of  the  variational 
prir  ).  Although  we  were  unable  to  determine  rigorous  upper  bounds  of  the  Hashin-Shtrikman 
and  n  varieties,  we  produce  some  arguments  justifying  the  application  of  the  approximate 

version  of  the  variational  procedure  in  connection  with  upper  bounds  for  the  linear  comparison 
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composite  to  obtain  "upper  estimates"  for  the  effective  properties  of  classes  of  nonlinear 
composites.  Finally,  the  question  of  attainability  of  the  bounds  is  discussed  briefly. 

In  the  development  of  the  bounds,  we  will  be  guided  by  the  physical  requirements  of  the 
deformation  theory  of  plasticity,  although  completely  analogous  results  will  also  hold  in  the  context 
of  nonlinear  creep.  Additionally,  we  could  follow  the  arguments  presented  by  DUVA  and 
Hutchinson  (1984)  to  make  use  of  the  effective  complementary-energy  function  of  a  given 
deformation-theory  composite  to  obtain  an  approximate  form  for  the  yield  function  of  a  flow- 
theory  composite  exhibiting  the  same  microstructure  (and  coinciding  plastic  response  under 
proportional  loading)  as  the  deformation-theory  solid. 

2.  Effective  Properties 

We  consider  a  heterogeneous  solid  occupying  a  domain  of  unit  volume  QeR^,  with 
boundary  dQ.  The  nonlinear  constitutive  behavior  of  the  solid  is  characterized  by  an  energy 
density  function,  H'(x,e),  depending  on  the  position  vector  x  and  the  strain  field  E(x),  in  such  a 
way  that  the  stress  field  0(x)  is  given  by 

(2.1)  Qi\)  =  d^w{\,e). 

Here  denotes  the  usual  differentiation  with  respect  to  e,  if  w  is  differentiable.  More  generally, 
can  be  given  the  interpretation  of  the  subdifferential  of  convex  analysis  (EKELAND  and  TEMAM, 
1974;  §  1.5). 

Following  Hill  (1963),  we  define  the  effective  constitutive  behavior  of  the  heterogeneous 
solid  in  terms  of  the  analogous  relation 

(2.2)  o  =  ^jVV'(e), 

where  o  denotes  the  mean  value  of  the  stress  field  over  Q,  and  W  refers  to  the  normalized  energy 
(recall  that  Q  has  unit  volume)  of  the  solid  when  subjected  to  the  uniform  boundary  condition 

(2.3)  u  =  ex,  xedQ, 
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where  u  is  the  displacement  field,  and  e  is  a  constant,  symmetric  tensor.  We  note  that  under  this 
boundary  condition,  the  mean  value  of  the  strain  over  Q  is  precisely  e. 

The  effective  energy  of  the  solid,  W ,  can  be  obtained  directly  in  terms  of  the  principle  of 
minimum  potential  energy  via  the  relation 

(2.4)  W(e)=  min  W(e), 

eeX(c) 

where 

(2.5)  W(e)  =  J  w(x,£(x))  dx 

Q 

is  the  pertinent  energy  functional,  and 

(2.6)  A'(e)  =  |e|  e  =  y[Vu  +  (Vu)^j  inQ,  and  u  =  ex  on 

is  the  set  of  kinematically  admissible  strain  fields. 

In  plasticity,  it  is  usually  more  convenient  to  express  the  constitutive  relation  of  the  solid  in 
terms  of  a  complementary-energy  density  function,  u  =  tv*,  where  >v’  is  the  Legendre  transform, 
or  polar  function  of  convex  analysis  (EKELAND  and  TEMAM,  1974;  §  1.4.1),  defined  by 

(2.7)  w*(x,a)  =  sup{ea- w(x,e)}. 

Thus,  the  constitutive  relation  of  the  heterogeneous  solid  in  its  local  form  can  alternatively  be 
written 

(2.8)  e(x)  =  ^„M(x,a). 

Correspondingly,  the  effective  constitutive  relation  of  the  solid  (in  its  global  form)  is  given  by 

(2.9)  z  =  d-U{o), 

where  0  corresponds  to  the  normalized  complementary  energy  that  would  exist  in  the 
heterogeneous  solid  when  subjected  to  a  uniform  constraint  on  dQ.,  such  that 

(2.10)  an  =  OT, 
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where  o  is  a  constant,  symmetric  tensor.  Thus,  the  effective  complementary  energy  of  the 
composite,  U ,  can  be  described  in  terms  of  the  principle  of  minimum  complementary  energy  via 
the  relation 

(2.11)  U(ja)=  mmU{G), 

oeS(a) 

where 

(2.12)  (/(a)  =  Ju(x.a(x))dx 

Q 

is  the  complementary  energy  functional,  and 

(2. 13)  S(a)  =  {a|  V  •  o  =  0  in  Q,  and  on  =  an  on  (?Q} 
is  the  set  of  statically  admissible  stress  fields. 

We  will  assume  that  w,  and  therefore  u,  are  convex  functions  of  e  and  O,  respectively,  and 
that  they  satisfy  certain  growth  conditions  to  be  specified  later.  Further,  by  convex  duality 
(EKELAND  and  TEMAM,  1974;  §1.4.2),  we  have  that 

(2.14)  w(x,e)  =  sup{eo-w*(x,a)}. 

Then,  the  effective  energy  functions  W  and  U  are  convex  (see  PONTE  CASTASeda  and  WILLIS, 
1988),  but  on  the  other  hand,  we  have  that,  in  general  (WILLIS,  1989), 

(2.15)  W(e)»s^{e- 0-^(0)}. 

The  reason  for  the  inequality  is  that  the  two  definitions  of  effective  properties  (in  terms  of  W  and 
0 )  correspond  to  different  boundary  conditions  (uniform  traction  versus  uniform  displacement). 
However,  for  a  composite  solid,  obtained  in  the  limit  as  the  size  of  the  typical  heterogeneity  in  the 
solid  becomes  small  compared  to  the  size  of  the  specimen,  the  two  definitions  agree,  and  equality 
holds  in  (2.15). 
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3.  The  New  Variational  Principles 

In  the  development  of  the  new  variational  principles,  we  will  consider  two  dual  formulations 
of  essentially  the  same  result,  depending  on  whether  we  start  from  the  minimum  potential  energy 
formulation  (2.4),  or  from  the  minimum  complementary  energy  formulation  (2.11).  In  both  cases, 
we  will  assume  that  the  solid  is  locally  isotropic,  so  that  we  can  write 

(3.1)  w(\,E)  =  0i\;€^,£j, 

where  x  R*  x  R*  R*  is  a  nonnegative  function,  that  is  convex  and  continuous  in  its  last 
two  arguments,  and  satisfies  the  condition  that  0(x;O,O)  =  0  for  all  x.  Here,  R*  is  the  set  of  the 
extended  nonnegative  reals  (R  is  the  set  of  extended  reals),  and 

(3.2)  e„=  jltrej  and  e,  = 

are  the  mean  and  effective  (in  their  plasticity  denotations)  strains,  respectively,  where  e  =  e  —  e„I 
is  the  deviator  strain  tensor.  We  note  that  the  above  conditions  on  0  are  consistent  with  the 
assumed  convexity  of  w,  and  that  the  form  (3.1)  is  not  the  most  general  form  for  the  energy 
function  of  a  nonlinear  isotropic  solid  (we  could  also  have  dependence  on  the  determinant  of  E), 
but  this  form  is  still  general  enough  to  cover  the  usual  plasticity  models  of  interest  here.  Also,  we 
emphasize  that  with  the  above  definitions  both  e„  and  e,  are  nonnegative. 

We  note  here  for  later  reference  that  under  the  above  hypotheses  on  w,  its  polar  function  can 
be  expressed  (see  EKELAND  and  Temam,  1974;  §  1.4.3  for  an  analogous  result,  and  also 
Appendix  I  for  a  definition  of  0*:Q  x  R*  x  R*  — >  R*)  in  terms  of 

(3.3)  vv‘(x,a)  =  <T,,3cr„  j, 

where 

(3.4)  (T„  =  ^|tral  and  cr.= 

are  respectively  the  mean  and  effective  stresses,  and  s  =  a  —  a„l  is  the  deviator  stress  tensor. 
Alternatively,  we  may  represent  w  =  w*  by 
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(3.5)  u(x,G)=y/(x;a^,aJ, 

where  ^  is  a  nonnegative  function,  convex  in  its  last  two  argument,  and  such  that 
Vr(x;0,0)  =  0  for  all  x.  These  properties  follow  from  simple  properties  of  the  polar  function  (see 
EKELAND  and  TEMAM,  1974;  §  1.4).  Finally,  we  have  that  w  =  u  can  also  be  represented  by 

(3.6)  u*(x,e)=  v^*|^x;|£„3e„j. 

3.1.  Minimum  potential  energy  formulation 

The  new  variational  principle  relies  on  a  quadratic  change  of  variables  u  =  h(e),  with 
h:R*  — >  R*  given  by  h(e)  =  e^.  We  can  then  define  a  function  /:fi  x  R*  x  R*  ^  R*  such  that 

(3.7)  /(x;u„u„)  =  0(x;r'(uJ,)i-‘(M„)). 

By  the  properties  of  0,  /is  a  nonnegative  function  with  same  dependence  on  x  as  0,  and  such  that 
/(x;0,0)  =  0.  We  note  further  that /is  continuous,  but  not  necessarily  convex  (in  its  last  two 
arguments).  In  fact,  because  we  are  interested  in  the  application  of  these  results  to  plasticity,  we 
will  usually  assume  that/is  concave,  or  at  least  that  it  has  concave  growth.  Thus,  for  example,  for 

w-fl 

an  incompressible  power-law  material,  w~e,"  (n  ^  1),  and  therefore  /~m/*,  which  is  a 
concave  function,  even  though  w,  itself,  is  convex. 

We  can  then  define  (see  Appendix  I,  and  VAN  TiEL,  1984;  §7.14)  the  concave  polar  function 
/.:OxR*xR*->R*  by 

(3.8)  /.(x;p.,p„)=  inf  {m.P,+m„P„-/(x;m.,m„)}, 
and,  therefore,  we  have  that 

(3.9)  f{x;u„uj<  inf  {u,p,  +  u„p„- fix-, 

P.-PmiO 

with  equality  if  /  is  concave. 

We  proceed  by  defining  the  functional  F  induced  by /  as  follows 

(3.10)  Fiu,,uJ  =  j  /(x;m,(x),m„(x))  dx, 

a 
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and  we  note  that  (see  EKELAND  and  Temam,  1974;  §  IX.2.1,  2)  that  the  (concave)  polar  of  the 
functional  F,  given  by 


(3.11)  FAp,,pJ=  inf  J  f[M,(x)/?,(x)  +  M„(x)p„(x)]tfe-F(u„u„) 


can  alternatively  be  represented  in  terms  of  /,  via 


(3.12)  F.(p,,p„)  =  j/.(x;p,(x),p„(x))  dr, 

Q 

with  an  analogous  representation  for  the  bipolar  function  F..  in  terms  of  /„ . 

Thus,  rewriting  the  minimum  energy  principle  (2.4)  in  the  form 

(3.13)  W{t)=  min  F(«.,u„), 

where  AT' is  the  set  of  admissible  fields  u^(x),  u„(x)  induced  by  K  in  (2.6),  we  arrive  by  means  of 
(3.9)  to  the  inequality 

(3.14)  Wi£)^W,(e), 
where 


(3.15) 


and  L  is  the  functional  given  by 


Q  Q 

We  emphasize  that  for  concave  f,  (e)  =  ll'(e) . 

Next,  we  simplify  the  expression  for  W+(e),  by  noting  that 
(3.17)  min  M  L{u„u^,p„pJ)=  inf  inf  L(m.,u„,p.,p„), 

u,.u,tK  P..P.20  P..P.20 

since  the  order  in  which  infima  are  evaluated  can  be  interchanged.  It  follows  that 


iv;(e)=  inf 

P.W.p.WiO 


w;(e)-  J  /.(x;p.(x).p„(x))  dxV 


(3.18) 
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where 


(3.19) 


^ «  ilL- 1 J +  w»(x)P*(x)]  dx  L 


3  2 

By  introducing  the  new  variables  =  — p,  and  K^=—p^,  and  expressing  and  in 

terms  of  e,  and  e„,  respectively,  we  are  able  to  rewrite  the  integrand  in  (3.19)  in  terms  of  the 
quadratic  energy  function 


(3.20) 


2  9 

w,(\,e)  =  <po{\;£„£„)  =  -fi,(\)e]+-}c^ix)£l. 


which  corresponds  to  a  linear-elastic,  isotropic,  heterogeneous  solid  with  nonnegative,  but 
otherwise  arbitrary,  shear  modulus  jUo(x),  and  bulk  modulus  x‘„(x).  With  this  notation,  the  new 

variational  principle  takes  the  final  form 
(3.21)  = 

where 


(3.22)  U;(e)  =  nun  f  w/x,e(x))  dx, 

eeX(c)J 

a 

and  where  V is  the  functional  specified  by  the  function  v(x;p„,  k^)  =  2 


(3.23) 


V'(P„  K'J  =  Jv(x;p^(x),  K-^(x))  dx. 


We  emphasize  once  again  that,  under  the  hypothesis  of  concave  /,  which  is  usually  valid  in 
plasticity,  the  new  variational  principle  and  the  classical  principle  of  minimum  potential  energy  are 
completely  equivalent,  and  we  have  equality  in  (3.14).  Thus,  under  these  circumstances,  the  new 
variational  principle  can  be  given  the  following  interpretation:  the  effective  behavior  of  the 
nonlinear  solid,  as  characterized  by  its  effective  energy  (2.4),  can  be  determined  alternatively  in 
terms  of  the  effective  energy  of  a  linear,  heterogeneous  comparison  solid,  with  variable  (within 
each  phase)  moduli,  p„(x)  and  v^(x),  whose  precise  variation  is  determined  by  the  new 


11 


variadcMial  principle  (3.21).  Clearly,  the  new  variational  principle  is  in  general  at  least  as  difficult  to 
evaluate  as  the  original  classical  variational  principle  for  the  nonlinear  solid.  The  main  advantage  of 
the  new  variational  principle  is  that  it  allows  useful  approximations  that  cannot  be  obtained  directly 
from  the  classical  variational  principle.  One  possible  approximation  in  (3.21)  is  to  replace  the 
innmum  over  the  set  of  (arbitrarily)  variable,  nonnegative  moduli  by  the  smaller  set  of  piecewise 
constant,  nonnegative  moduli.  We  will  consider  the  application  of  this  approximation  in  the 
determination  of  bounds  for  the  effective  properties  of  nonlinear  composites  in  the  following 
section.  This  type  of  approximation  was  introduced  directly  by  PONTE  CaSTANTEDA  (1991a).  In 
particular,  it  was  found  in  that  work  that  when  this  approximation  was  used  in  conjunction  with  the 
well-known  Hashin-Shtrikman  bounds  for  linear  composites,  nonlinear  bounds  resulted  that  were 
in  agreement  (see  also  WiLLIS,  1991)  with  earlier  bounds  determined  earlier  by  PONTE 
Castaneda  and  Willis  (1988)  using  a  nonlinear-extension  of  the  Hashin-Shtrikman  variational 
principles  due  to  WILLIS  (1983)  and  Talbot  and  WiLLlS  (1985).  We  note  that  the  new  variational 
principle  is  restricted  [for  equality  to  hold  in  (3.14)]  thus  far  to  isotropic  phases,  although  it  can  be 
used  in  conjunction  with  anisotropic  microgeometries  leading  to  anisotropic  effective  properties. 
On  the  other  hand,  the  Talbot-Willis  variational  principle  can  be  used  to  obtain  Hashin-Shtrikman 
bounds  for  composites  with  anisotropic  phases,  although  some  level  of  approximation  is  required 
for  computational  convenience  (see  DENDIEVEL,  BONNET  and  WILLIS,  1991). 

Finally,  we  remark  that  statement  (3.9)  can  be  rewritten  in  the  form 

(3.24)  w(x;e)<  inf  {w^(x,e)- v(x;/i,,fc-j}, 

(with  equality  if  w  is  such  that/ is  concave),  where  v(x;/i„,  tc^)  takes  the  corresponding  form 

(3.25)  v(x;/t„,K-„)  =  sup{H'(x,e)-  w„(x,e)}. 

With  this  representation  of  the  local  energy  functions,  the  new  variu.ional  principle  is  seen  to  be 
nothing  more  that  the  global  version  of  the  relation  (3.24).  Although  we  have  already  remarked  that 
the  new  variational  principle  (3.21)  is  only  exact  for  isotropic  phases,  provided  that/is  concave, 
the  form  (3.25)  for  v(x;/i„,  v^)  suggests  that  the  new  principle  could  al.^o  be  used  for  composites 
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with  anisotropic  phases.  Indeed  this  is  the  case,  but  in  this  event,  equality  could  not  hold  in  (3.24), 
and  the  new  variational  principle  in  the  given  form,  would  yield  only  an  upper  bound 
approximation  to  the  effective  behavior  of  the  nonlinear  heterogeneous  solid  W,  that  would  be 
expected  to  get  progressively  weaker  with  the  degree  of  anisotropy  in  the  constituent  phases. 
Improved  bounds  could  be  obtained  by  the  introduction  of  anisotropic  comparison  materials.  This 
is  outside  the  scope  of  this  paper,  and  will  be  pursued  elsewhere. 

3.2.  Minimum  complementary  energy  formulation 

Given  the  function  yr  defined  by  (3.5),  we  make  use  of  the  same  change  of  variables  h  as  in 
the  previous  subsection  to  define  a  function  x  R*  x  R*  — >  R*  such  that 

(3.26)  ^(x;v.,v„)= 

By  the  properties  of  y/,  enunciated  at  the  beginning  of  this  section,  g  is  a  nonnegative  function  with 
g(x;0,0)  =  0.  Further,  g  is  continuous,  but  not  necessarily  convex,  in  its  last  two  arguments. 
However,  it  is  shown  in  Appendix  II  that  if/ is  concave,  then  g  is  convex.  As  mentioned  earlier, 
this  is  the  usual  situation  in  plasticity.  For  instance,  for  the  same  example  of  an  incompressible 

ii-fl 

power-law  material  given  in  the  previous  subsection,  u  ~  cr"*‘  (n  >  1),  and  g  ~  ,  which  are 

both  convex  functions. 

We  can  then  define  (see  Appendix  I)  the  convex  polar  function  g*:Q  x  R*  x  R*  — >  R*  via 

(3.27)  sup  +  v„i7„ -g(x;v„v„)}, 

so  that 

(3.28)  g(x;v„v„)>  sup  K<7,  +  v„i7„-g‘(x;(7„q„)}, 

with  equality  if  g  is  convex. 

Introducing  the  functional 

(3.29)  G(v.,v„)  =  Jg(x;v,(x),v„(x))dx, 

a 
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and  following  a  procedure  analogous  to  that  of  the  previous  subsection,  we  arrive  at  the  following 
inequality  involving  the  effective  complementary  energy  of  the  nonlinear  heterogeneous  solid  as 
characterized  by  (2.11),  namely, 

(3.30)  U(a)>U_(a), 


where 


(3.31) 


U.(o)=  min  sup  A/(v.,v„, 


M  is  the  "saddle"  functional  given  by 


(3.32)  Af(v.,v.,<7.,<7^)  =  J[v,(x)<7,(x)  +  v„(x)<7„(x)]dx-Jg‘(x;<7.(x),<7.(x))dx, 

a  a 

and  S'  is  the  set  induced  by  the  above  change  of  variables  on  the  set  of  admissible  stresses  S.  We 
emphasize  that  for  concave /(and  therefore  convex  g),  U_(a)  =  U  (o) . 

In  order  to  simplify  the  above  expression  for  U_,  we  are  required  to  interchange  the  order  of 

the  minimum  and  the  supremum.  This  is  allowed  by  Prop.  VI.2.3  of  EKELAND  and  Temam 
(1974),  because  M is  affine  and  therefore  convex  in  v,,  v^,  and  concave  (because  g*  is  convex)  in 

<7,,  q^,  and  lim  =  ©o  for  q^°^  >  0  fixed.  Thus,  we  conclude  that 


(3.33) 


min  sup  M{v„v„,q,,q„)=  sup  inf  Af(v,,v^,<7,,<7^). 


This  leads  to  the  following  restatement  of  (3.31) 

(3.34)  0.(a)  =  sup  W,(o)  -  V(/z„,  kJ}, 


where 


(3.35) 


f7,(o)=  min  fM„(x,a(x))dx 

0€5(3)J 

a 


is  the  effective  complementary  energy  of  the  comparison  composite,  with  =  w*  [see  (3.20)],  and 


/  1  1  "1 

where  the  functional  V  is  given  by  relation  (3.23),  with  v(x;^„,  k^)  =  g*  x;- — ,- — 

2k J 


.  This  last 


expression  for  v  can  be  expressed  directly  in  terms  of  u  and  via 
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(3.36)  v(x;/i„,  K-J  =  sup{m„(x,o)  -  M(x,a)} . 

In  terms  of  this  representation,  it  is  easy  to  verify  that  the  v  functions  given  by  (3.25)  and  (3.36) 
are  in  fact  identical  (and  hence  the  choice  of  the  same  names).  Simply  note  that  from  (3.36),  and 
the  fact  that  the  order  of  suprema  can  be  interchanged,  it  follows  that 

v(x;/t^,  K-J  =  sup|sup{e  •  a-  w„(x,e)}  -  M(x,a)| 

37)  =sup{sup{e  o-«(x,o)}-w^(x,e)J 

=  sup{w(x,e)- w„(x,e)}. 


which  is  in  agreement  with  (3.25).  This  result  can  also  be  obtained  by  means  of  result  (II.7)  in 


Appendix  n,  which  shows  that  f.{p,,pj)  =  -g* 


'  1  _9 


.  Also,  note  the  following  counterpart 


to  (3.24) 

(3.38)  u(\,a)>  sup  {M„(x,CT)-v(x;/t<„vj}, 

which  is  the  local  version  of  (3.34). 

Finally,  we  study  the  relation  between  the  two  variational  principles  and  (/_,  as  given  by 
(3.21)  and  (3.34).  We  begin  by  evaluating  the  polar  of  yielding 


(3.39) 


w;(a)  =  sup{ea-W^(e)} 

=  supjeo-  inf  {MC(e)  + V'(/i„,K-j}} 
=  r  ,Jsup{e  a-M>(e)}-V(/i„,x-j| 

/<,(»), >f.(x)20l  cl  j  J 

=  sup  {w;(a)-V(/i„,K-j}, 


where  once  again  we  have  used  the  fact  that  the  order  of  the  suprema  can  be  interchanged.  Next, 
by  specializing  the  dual  version  of  result  (2.15)  to  the  linear  comparison  composite  (/.«., 
^  (7,),  we  conclude  that 
(3.40) 
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For  a  composite  (see  the  discussion  following  (2.15)),  =  0^ ,  and  therefore  we  have  W’  =  U_ 

in  place  of  (3.40),  which  demonstrates  that  there  is  no  duality  gap  between  the  minimum  potential 
and  minimum  complementary  energy  versions  of  the  new  variational  principle.  Further,  if/ is 

concave,  then  we  have  equality  in  (3.14)  and  (3.30),  simultaneously,  and  therefore  we  also  have 
equality  between  W'  and  U ,  which  is  in  agreement  (as  it  should)  with  the  discussion  following 

(2.15).  Hence,  in  the  application  of  the  new  variational  principle,  one  may  utilize  either  one  of  the 
two  versions  of  the  variational  principle.  The  choice  of  the  version  will  usually  be  dictated  by 
convenience. 


4.  Application  to  Bounds 

In  this  section,  we  are  concerned  with  the  determination  of  bounds  for  the  effective  properties 
of  composite  materials  in  the  context  of  plasticity.  The  approach  will  be  to  make  use  of  well-known 
bounds  for  linear-elastic  composites  in  the  new  variational  principles  to  induce  corresponding 
bounds  on  the  effective  properties  of  nonlinear  plastic  composites.  Although  the  procedure  can  be 
applied  more  generally  than  we  will  consider  in  this  section,  we  provide  a  sampling  of  the  possible 
applications  of  the  method  by  including  a  new  derivation  of  the  classical  Voigt/Reuss  bounds  for 
generally  anisotropic  composites,  new  forms  for  the  Hashin-Shtrikman  bounds  for  isotropic  and 
transversely  isotropic  composites,  and  completely  new  higher-order  bounds  of  the  Beran  type  also 
for  isotropic  composites.  We  will  be  interested  only  in  materials  with  energy  functions  satisfying 
the  concavity  of/hypothesis,  which  is  quite  adequate  in  plasticity.  In  this  event,  we  have  complete 
equivalence  between  the  new  and  classical  variational  principles,  and  it  is  then  clear  from  either 
(3.21)  or  (3.34)  that  bounds  on  the  effective  energy-density  functions  of  the  linear  comparison 
composites  can  be  utilized,  at  least  in  principle,  to  yield  corresponding  bounds  on  the  effective 
energy-density  functions  of  nonlinear  composites. 

In  particular,  we  will  study  heterogeneous  solids  with  n  distinct  phases,  each  of  which  is 
homogeneous  and  isotropic.  The  volume  fractions  of  each  phase  will  be  assumed  to  be  given. 
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According  to  relation  (3.1),  we  can  characterize  the  energy-density  functions  of  the  different 
phases  in  terms  of  functions  x  R*  ->  R*  (r  =  1, ....  n),  such  that  the  local  energy-density 

function  of  the  composite  is  given  by 

n 

(4.1)  0(x;e„e„)= 

r=l 

where  ;t^'^(x)  is  the  characteristic  function  of  phase  r  (this  function  vanishes,  unless  x  is  in  phase 
r,  in  which  case  it  equals  unity).  However,  because  of  the  usual  representation  in  plasticity  of  the 
total  strain  as  the  sum  of  the  (linear)  elastic  strain  and  the  (nonlinear)  plastic  strain,  we  find  it  more 
convenient  to  represent  the  constitutive  relation  for  the  composite  in  terms  of  the  corresponding 
complementary  energy-density  functions  of  each  phase,  with  a  representation  for  the 
complementary  energy -density  function  of  the  composite  y/ completely  analogous  to  (4.1).  We  will 
restrict  our  attention  to  the  class  of  complementary-energy  functions 

(4.2)  + 

where  the  <p^'^  are  convex,  nonnegative  functions  (and  such  that  the  corresponding/ functions  are 
concave),  and  are  nonnegative  given  constants.  The  form  (4.2)  is  quite  appropriate  in 
plasticity  since  the  plastic  strain  is  usually  taken  to  be  traceless  (incompressible),  so  that  the 
compressible  part  of  the  constitutive  relation  is  due  only  to  the  elastic  strain,  and  is  hence  linear 
(corresponding  to  a  quadratic  potential).  Thus,  is  the  (elastic)  bulk  modulus  of  phase  r,  and 
<p^'^  characterizes  the  total  shear  behavior  of  phase  r  (both  elastic  and  plastic  contributions).  For 
instance,  for  a  linear-elastic  material, 

(4.3) 

where  ^  0  is  the  elastic  shear  modulus  of  the  material.  More  generally,  a  plastic  potential  will 
usually  be  added  to  the  potential  in  (4.3)  to  account  for  the  plastic  strains.  In  our  development  in 
the  following  sections,  we  will  not  specialize  our  consideration  to  any  specific  form  for  the  plastic 
potential,  and  our  results  will  hold  for  arbitrary  energy-density  functions  of  the  form  (4.2),  with 
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the  given  hypotheses  on  and  k^'^.  We  further  note  that  the  volume  fractions  of  each  phase  arc 
given  by 

(4.4)  dx 

a 

W 

and  are  such  that  =  1 . 

r=\ 

In  accord  with  the  above-mentioned  representation  for  the  local  complementary  energy- 
density  function  u{x,<s)=  we  will  make  use  of  the  complementary  energy 

formulation  (3.34)  of  the  new  variational  principle,  although  equivalent  results  would  be  obtained 
by  means  of  the  dual  potential  energy  formulation  (3.21).  Thus,  we  will  consider  first  in  §  4.1 
lower  bounds  for  the  effective  energy  function  U ,  and  then  continue  in  §4.2  with  the  harder 
problem  of  determining  upper  bounds  for  U.  The  relationship  of  the  bounds  to  actual 
microstructures,  and  concomitant  implications  for  optimality  of  the  bounds  will  be  addressed 
briefly  in  §  4.3. 

4. 1 .  Lower  bounds 

As  indicated  in  §  3,  lower  bounds  for  the  effective  complementary  energy  of  the  nonlinear 
composite  0  can  be  obtained  directly  from  (3.34)  by  restricting  the  class  of  completely  arbitrary 
(nonnegative)  comparison  moduli  to  the  smaller  class  of  piecewise  constant  (nonnegative)  moduli. 
Thus,  in  the  case  of  composites  characterized  by  (4.2),  we  consider  a  comparison  composite  with 
shear  moduli  given  by 

(4.5)  = 

rsl 

where  the  are  nonnegative  constants.  We  note  that  we  have  selected  the  same  microstructural 
distribution  for  the  comparison  composite  as  for  the  nonlinear  composite.  Also,  we  remark  that 
because  the  "compressive"  terms  in  the  functions  are  quadratic,  the  supremum  problem  over 
the  comparison  bulk  moduli  in  (3.34)  is  solved  exactly  by 
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(4.6)  = 

r=I 

where  K^p  -  k^'K 

It  follows  from  this  discussion  that 

(4.7)  U(c)  >  ^sup^|(7.(a)  - 

where  U„  corresponds  to  the  effective  complementary  energy  of  the  linear  comparison  composite 

with  shear  and  bulk  moduli  characterized  by  (4.5)  and  (4.6),  respectively,  and  where  the  functions 
take  the  simplified  form 

(4.8)  v^'^(uP)= ~  ?>^'^(-y)j- 
We  note  for  later  reference  that 


(4.9) 


Thus,  any  lower  bound  for  0^  can  be  utilized  in  (4.7)  to  induce  a  lower  bound  on  U.  We 


note  that  for  nonlinear  composites,  just  as  for  linear  composites,  the  bounds  on  the  effective  energy 
will  depend  on  the  prescribed  microstructural  information.  For  example,  the  Hashin-Shtrikman 
bounds  for  isotropic  composites  will  be  tighter  than  the  Voigt/Reuss  bounds  for  generally 
anisotropic  composites.  However,  it  is  important  to  note  that  the  characterization  of  microstructure 
is  different  for  linear  and  nonlinear  materials.  For  example,  a  periodic  linear  composite  with  unit 
cell  possessing  icosahedral  symmetry  leads  to  isotropic  effective  behavior  for  the  composite.  The 
same  is  not  true  in  general  of  a  periodic  nonlinear  composite  with  an  icosahedral  unit  cell.  This  is 
because  linear  behaviors  can  be  characterized  by  (finite  dimensional)  tensors,  whereas  nonlinear 
behaviors  require  in  general  characterization  in  terms  of  (infinite  dimensional)  functions. 
Notwithstanding  this  difficulty,  we  observe  that  the  class  of  isotropic  linear  microstructures 
includes  the  class  of  isotropic  nonlinear  microstructures.  Hence,  nonlinear  bounds  based  on  linear 


isotropic  microstructures  will  also  hold  for  nonlinear  isotropic  microstructuies.  Analogous  results 
would  hold  for  any  other  material  symmetries  for  the  composite  (e.g.  transversely  isotropic).  Thus, 
a  nonlinear  lower  bound  for  nonlinear  composites  with  a  given  symmetry  results  by  replacing 
in  (4.7)  by  a  lower  bound  for  with  the  same  symmetry. 

4.1.1.  Voigt  bound.  We  begin  by  considering  the  case  of  a  generally  anisotropic  composite 
with  phase  potentials  (4.2)  in  given  volume  fractions  Thus,  in  this  case,  we  require  a  lower 
bound  for  the  effective  complementary-energy  function  of  a  linear  comparison  composite  with 
shear  and  bulk  moduli  characterized  by  (4.5)  and  (4.6),  in  volume  fractions  c^'\  Such  a  bound  is 
provided  by  the  so-called  Voigt  estimate 


(4.10) 


1  ^2  ,  1  •s:2 


(/r'(o)= 


cr;+- 


w  n 

where  and  =  This  result  can  be  readily  obtained  from  the 


r=l 


principle  of  minimum  complementary  energy  (Hill,  1952;  PAUL,  1960).  Thus,  use  of  this  lower 
bound  for  the  linear  comparison  composite  into  (4.7)  leads  to  the  following  lower  bound  for  the 
effective  energy  of  the  nonlinear  composite 


(4.11) 


1  j=2 

■<y_i 


r=l 


2k' 


(V) 


where  The  above  bound  involves  2n  optimizations  (two  for  each  phase),  but  it 

r«l 


can  be  greatly  simplified  by  means  of  the  identity  of  Appendix  Ill 


20 


where  the  infimum  is  over  the  set  of  variables  (r  =  1,...,  n),  which  are  subject  to  a  zero- 

m 

average  constraint  a  =  =  0  .  With  this  identity,  the  expression  of  the  bound  in  (4.11) 

r=l 

leads  to 

(4.13)  U^''\a)=  sup<  inf 

a)*'* 

S)=0 

Next,  recalling  that  is  concave  in  (since  g*  is  convex),  and  noting  that  the  expression 

inside  the  square  brackets  is  convex  in  o)^'\  the  Saddle  Point  Theorem  (EKELAND  and  TEMAM, 
1974;  §  VI.4.2)  allows  the  interchange  of  the  order  in  which  the  suprema  and  infima  are  evaluated 
to  arrive  at 

(4.14)  =  lie"’  1+2^^"-’ 

3=0  ^ 

which  reduces  finally  to 

(4.15)  =  inf 

5J=0 

where  we  have  made  use  of  (4.9).  To  the  knowledge  of  the  author,  this  result  is  new. 

This  lower  bound  for  U,  obtained  from  the  linear  Voigt  bound  via  the  new  variational 
principle,  can  be  shown  to  be  equivalent  to  the  well-known  bound,  obtained  directly  from  the 
principle  of  minimum  complementary  energy  (for  a  nonlinear  composite)  by  assuming  a  constant 
stress  field  over  the  composite, 

(4.16)  (/(g)2fXc‘'’(4>''')''|  {S.)  +  :^Sl 

v--=i  y 

For  this  reason,  we  refer  to  the  lower  bound  (4. 1 .5)  as  the  nonlinear  Voigt  bound.  We  note  that  the 
standard  form  of  the  bound  (4.16)  involves  n  +  1  optimizations  (Legendre  transforms),  whereas 
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the  new  forai  (4.15)  involves  n  optimizations  with  a  simple  constraint  For  the  case  of  two  phases, 
the  new  fram  of  the  Voigt  bound  reduces  further  to 

(4.17)  U‘''>(o)  =  inf{c“V'''(|l-c“’al|^.)+c«V™(|l+c“><S)l^,)}+-^5^ 


which  involves  only  one  optimization. 

4.1.2.  Hashin-Shtrikman  lower  bounds.  Next,  we  consider  nonlinear  composites  with 
phases  having  complementary  energy-density  functions  (4.2),  in  prescribed  volume  fractions 
distributed  in  such  a  way  that  the  overall  composite  is  isotropic.  In  this  case,  the  linear  comparison 
composite  to  be  used  in  conjunction  with  (4.7)  is  taken  to  be  isotropic  with  shear  and  bulk  moduli, 
given  once  again  by  (4.5)  and  (4.6),  respectively,  in  volume  fractions  and  with 

The  relevant  isotropic,  linear  lower  bound,  due  to  HASHIN  and  SHTRIKMAN  (1963),  is  provided 
by  the  relation 


(4.18) 


1  _2  .  1-2 


where 
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where  we  have  replaced  by  in  the  expression  (4.19)2  for  Note  that  this  expression 

still  depends  on  fil'K  Similarly,  we  replace  by  =  max|K‘^''^|  in  expression  (4.19)i  for 

„(«S+) 

r*o 

In  general,  the  expression  (4.20)  for  the  nonlinear  isotropic  lower  bound  is  quite 
complicated.  Simplification  can  be  achieved  by  rewriting  the  above  expression  for  in  the 

form 


(4.21) 


which  in  turn  can  be  rewritten  in  the  form 


(4.22) 


by  means  of  the  identity  of  Appendix  III.  Similarly,  we  may  write 


(4.23) 


fsQ  ^  ^ 


ft 

where  the  opdmizadon  variables  satisfy  the  constraint  y  =  =  0.  Therefore,  we  have 


t/<“->(a)=sup|  inf 

y">.f=0  ~  '' 


(4.24) 


By  the  same  arguments  following  (4.13),  allowing  the  interchange  of  the  suprema  and  infima,  we 


amveat 
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(4.25) 


^<®-)(o)=  inf 


■■■+ /'’f  o-= 


Making  use  of  the  expression  for  this  can  in  turn  be  rewritten 


f 

r  /•  f 

_  _ 

i7‘"^->(a)  =  mim 

1 

inf  1 

sup- 

Uy'lo 

S‘" 

.Md  H’o 

y’KfmO 

L 

(4.26) 


where  we  have  used  the  fact  that  -4;r  =  ^ 


is  a  monotonically  decreasing  function 


of  This  result  finally  leads  to 
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with 
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(4.29)  s  = 
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Note  that  s  depends  on  and  therefore  the  supremum  in  (4.28)  is  not  trivially  evaluated.  Thus, 

the  new  expression  for  the  i.  .niinear  Hashin-Shtrikman  lower  bound  involves  an  2/i-diniensional 
constrained  optimization  as  given  by  (4.27),  together  with  the  optimization  problem  for  Also, 

we  note  that  there  are  n  "branches"  to  the  solution  (one  for  each  phase),  and  the  minimum  over  all 
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the  branches  yields  the  desired  lower  bound.  Additionally,  we  remark  that  the  expression  for  the 
bound  simplifies  to 


if  at  least  one  phase  is  incompressible  (so  that  —>««).  If  all  phases  are  incompressible,  the 
expression  simplifies  further  to 


(4.32) 


^(«s-)(0)  _  1  inf  J  f  f’ 


where  the  are  the  same  as  above,  and 


(4.33)  = 

If  there  are  only  two  phases,  the  constraints  can  be  accounted  for  by  letting  and 

(with  similar  forms  for  and  /^’).  Then,  when  both  phases  are  incompressible, 
we  arrive  at  the  simple  expression 


25 


This  expression  generalizes  results  given  by  PONTE  Casta^DA  (1991b)  for  a  dual-phase 
incon^nessible  composite  with  a  linear  and  a  nonlinear  phase.  Another  simple  result  is  provided  by 
the  porous  material,  when  one  phase,  say  phase  2,  is  void  (the  other  is  incompressible).  Then,  the 
lower  bound  (4.30)  reduces  to  (since  - 1) 

(4.35)  (7(™-)(a)  =  c<*V^'f47 

first  given  by  PONTE  Casta^DA  (1991a).  WILLIS  (1991)  has  also  shown  explicitly  that  this 
particular  result  can  be  obtained  from  the  Talbot-Willis  variational  principle  by  proper  choice  of  the 
homogeneous  comparison  material.  A  weaker  version  of  this  result,  using  the  Talbot-Willis 
variational  principle,  was  provided  by  PONTE  Castaneda  and  WiLLIS  (1988).  'ITiere  are  other 
simple  cases  in  the  context  of  isotropic  nonlinear  composites,  but  in  the  interest  of  time,  they  will 
be  considered  elsewhere. 

Finally,  we  remark  that  the  new  variational  principle  can  also  be  used  in  connection  with 
bounds  for  linear  anisotropic  composites  with  special  symmetries  to  obtain  corresponding  bounds 
for  nonlinear  anisotropic  composites  with  given  symmetries  (the  first  results  given  in  this  section 
hold  for  generally  anisotropic  composites).  For  instance,  the  case  of  fibre-reinforced  materials  is 
considered  by  PONTE  CastaNeda  and  De  BOTTON  (1991).  These  authors  make  use  of  the 
corresponding  linear  bounds  due  to  Hill  (1964),  and  Hashin  (1965)  (see  also  WALPOLE,  1969). 
The  results  are  more  complicated  than  the  results  for  the  isotropic  composite  (since  they  have 
transversely  isotropic  symmetry),  but  have  essentially  the  same  form.  For  illustrative  purposes,  we 
include  the  final  result  for  a  composite  with  incompressible  phases 

(4.36)  ^^""->(o)  =  mm-  inf  + 


with 
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are  the  transversely  isotropic  invariants  of  the  mean  stress  o  implied  by  the  fiber  orientation  vector 
n.  Results  in  a  different  form  for  this  case  have  also  been  given  by  TALBOT  and  WILLIS  (1991), 
using  the  Talbot-Willis  variational  method. 


4.2.3.  Higher-order  lower  bounds.  So  far,  we  have  given  examples  of  bounds  that  are  either 
well-known,  or  can  be  obtained,  at  least  in  principle,  by  other  procedures.  Thus,  the  Voigt 
nonlinear  bound  can  be  obtained  more  directly  from  the  classical  principle  of  minimum 
complementary  energy,  and  the  nonlinear  Hashin-Shtrikman  lower  bounds  can  also  be  obtained  by 
means  of  the  Talbot-Willis  nonlinear  extension  of  the  Hashin-Shtrikman  variational  principle. 
However,  the  simple  forms  (4.15)  and  (4.17)  for  the  Voigt  bound,  and  (4.27)  to  (4.32)  for  the 
Hashin-Shtrikman  bounds  are  new.  In  this  subsection,  we  discuss  the  determination  of  nonlinear 
bounds  that  have  not  been  possible  to  determine  by  any  other  method  hitherto. 

Although  many  different  higher-order  bounds  are  available  for  linear-elastic  composites,  we 
will  restrict  our  consideration  to  the  third-order  bounds  of  BERAN  and  MOLYNEUX  (1965),  and 
McCoy  (1970).  The  reason  is  that  while  higher-order  bounds  can  have  quite  complex  forms 
involving  many  additional  geometric  parameters  (other  than  the  volume  fractions),  MILTON  (1982) 
has  given  a  particularly  simple  form  for  the  third-order  bounds  of  dual-phase  isotropic  composites, 
which  involve  only  two  additional  geometric  parameters.  Milton's  form  for  the  bounds  and 
may  be  expressed  by 


(4.38) 


-1 


1-1 
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which  is  identical  in  form  to  (4.21),  but  with  jj}*^  =  0/6,  where  6  is  given  by  a  ratho*  complicated 
expression  [see  (3.29)  in  MILTON  (1982)]  involving  and  as  well  as  the  geometric 
parameters  and  and  by 


(4.39) 


1 


4  * 

with  =  .  Here,  the  geometric  parameters  and  lie  in  the  interval  [0, 1],  and 

^  r-I 


Substitution  of  these  expressions  for  the  comparison  moduli  into  (4.7)  yields  an  isotropic 
lower  bound  for  the  effective  energy  of  the  nonlinear  composite  involving  the  phase 
potentials  and  of  the  two  phases,  the  corresponding  volume  fractions  and  and 
the  two  additional  geometric  parameters  rf^  and  However,  due  to  the  complicated  form  of 
in  (4.38),  no  further  simplification  of  the  nonlinear  bound  is  available.  Fortunately,  for  the 
case  of  two  incompressible  phases,  the  expression  for  simplifies  to 


(4.40) 


r*l 

0) 


which  depends  only  on  the  geometric  parameter  rf'.  Then,  use  of  a  procedure  very  similar  to  the 
one  following  (4.21)  for  the  Hashin-Shtrikman  lower  bound  leads  to  the  nonlinear  third-order 
lower  bound 


6r'“-»(a)  =  inNc''V 

«».r 


^(1  -  o,  V— 


(4.41) 


IV  3  , 


where  =  1  -  We  remark  that  the  the  corresponding  nonlinear  Hashin-Shtrikman  lower 
bound  (4.34)  follows  immediately  from  this  it  by  choosing  ier  -  0,  or  77^'’  =1, 


whichever  yields  the  lowest  value  (note  that  the  infimum  problem  over  y  disappears  in  either  case). 
This  is  completely  analogous  to  the  corresponding  result  for  linear  two-phase  incompressible 
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composites  (MILTON,  1982).  Finally,  we  note  that  a  superficially  simpler  result  is  obtained  by 
making  use  of  the  potential  energy  formulation  instead  of  the  complementary  energy  formulation, 
as  we  have  done  above,  since  the  infimum  over  y  does  not  appear  in  the  corresponding  form  for 
However,  such  an  expression  will  in  general  be  harder  to  evaluate,  because  in  plasticity  the 
form  of  the  energy  function  of  the  phases  is  generally  more  complicated  than  the  corresponding 
form  for  the  complementary  energy  functions,  used  in  (4.41). 


4.2.  Upper  bounds  and  estimates 

It  is  clear  that  the  approximate  version  of  th«  new  variational  principle  (4.7)  will  not  be 
helpful  in  the  determination  of  upper  bounds  for  U .  Therefore,  we  must  resort  to  the  exact  version 
of  the  new  variational  principle  (3.34),  which  makes  the  determination  of  upper  bounds  for  U  an 
intrinsically  harder  problem  than  the  corresponding  lower  bound  problem.  In  this  subsection,  we 
will  attempt  to  go  as  far  as  we  can  with  the  problem  of  rigorous  upper  bounds,  and  then  revert  to 
the  simpler  problem  of  finding  estimates  for  the  upper  bound,  or  "upper  estimates".  This  is 
accomplished  essentially  by  making  use  of  the  approximation  (4.7)  in  spite  of  the  difficulty  just 
mentioned.  In  §  4.3,  we  will  provide  some  justificauon  for  making  this  approximation,  and  give  an 
interpretation  for  the  "upper  estimates"  for  0 . 

4.2.1.  Reuss  bound.  We  begin  by  considering  the  generally  anisotropic  composite  made  up 
of  phase  energy  functions  (4.2)  in  prescribed  volume  fractions  without  any  further  restriction 
of  the  microstructure.  Then,  the  effective  complementary  energy  function  of  the  linear  comparison 
composite  with  arbitrarily  inhomogeneous  moduli  /t„(x)  and  v„(x)  can  be  shown  (by  means 

of  the  principle  of  minimum  complementary  energy)  to  be  bounded  above  by 


(4.42) 
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(4.43) 


and 
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These  expressions  are  the  continuous  versions  of  the  standard  (discrete)  Reuss  estimates. 

Correspondingly,  the  effective  complementary  energy  function  of  the  nonlinear  composite  U 
is  bounded  above  by 


(4.44) 


2jc' 


where 


(4.45) 


v(x;// J  =  supj^/  -  «>(x;j)|. 


and  where  = 


■  fXr) 
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,  since  the  hydrostatic  strain  of  the  nonlinear  composite  is  taken  to  be 


linear.  Additionally,  we  have  that 


(4.46) 


which  allows  further  simplification  of  (4.44),  since 


(7<*\o)=  sup  \ 
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and  therefore 
(4.48) 


This  last  result  is  the  Reuss  estimate  for  the  nonlinear  composite,  which  could  have  been  obtained 
from  principle  of  minimum  complementary  energy,  applied  directly  to  the  nonlinear  composite. 
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The  present,  less  direct,  derivation  of  this  result  via  the  new  variational  principle  serves  to 
demonstrate  that  the  new  method  can  at  least  recover  the  classical  bounds  exactly.  This  derivation 
serves  further  to  illustrate  how  the  new  variational  principle  may  be  used  in  other  contexts  for 
which  upper  bounds  may  not  be  available  via  the  classical  principles. 

Interestingly,  we  remark  that  the  above  result  for  the  nonlinear  Reuss  bound  (4.48)  could 
also  have  been  obtained  by  means  of  the  approximation  (4.7),  following  a  procedure  analogous  to 
that  used  in  the  derivation  of  the  nonlinear  Voigt  bound  (4.15).  This  suggests  that  the 
approximation  introduced  in  (4.7)  is  exact  in  this  case.  We  will  give  some  insight  in  the  reasons 
behind  this  phenomenon  in  §  4.3. 

4.2.2.  Hashin-Shtrikman  upper  estimates.  We  consider  next  the  upper  bound  problem  for 
the  class  of  nonlinear  composites  with  phase  potentials  (4.2)  in  fixed  volume  fractions  (4.4), 
distributed  in  such  a  way  that  the  effective  behavior  of  the  composite  is  isotropic.  In  this  case, 
however,  direct  application  of  the  exact  version  of  the  new  variational  principle  (3.34)  leads  to  a 
much  harder  problem  than  the  one  solved  above  for  the  Voigt  bound  for  anisotropic  geometries. 
For  this  reason,  we  will  not  pursue  this  approach  here,  and  resort  instead  to  the  use  of  the 
approximation  introduced  in  (4.7).  Then,  upper  bounds  for  the  effective  energy  0^  of  the  class  of 

isotropic  linear  comparison  composites  with  prescribed  volume  fractions  lead  to  estimates  for  the 
upper  bound  of  the  class  of  isotropic  nonlinear  composites  with  prescribed  volume  fractions.  We 
will  refer  to  these  estimates  as  "upper  estimates"  since  they  may  be  close  to  the  upper  bound.  (In 
fact,  the  corresponding  computation  for  the  Reuss  upper  oound  suggests  that  in  some  cases  these 
estimates  may  yield  rigorous  bounds).  Further,  in  practical  applications,  we  are  rarely  interested  in 
precise  bounds;  instead,  estimates  are  often  preferred. 

The  computation  of  the  "upper  estimates"  follows  a  procedure  completely  analogous  to  the 
computation  of  the  lower  bounds.  Thus,  the  effective  energy  function  of  the  linear  comparison 
composite  0^  is  bounded  above  by 
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where  are  given  by  expressions  similar  to  (4.19),  except  that  and 


must! 


be  replaced  by  and  =  nun{Ki''’}.  Then,  use 


of  this  upper  bound  for  into  (4.7)  leads  to  the  "upper  estimate"  for  U .  The  final  result 
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where  is  as  given  by  (4.28)  and  (4.29),  but  with  replaced  by  x*^  \  Note  that,  apart  from 
these  modifications,  the  expression  (4.50)  for  differs  from  the  expression  (4.27)  for 


only  in  that  the  outermost  minimum  over  all  phases  has  been  replaced  by  a  maximum. 

When  all  the  phases  are  incompressible  (so  that  the  above  result  simplifies 

further  to 


(4.51) 
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where  and  are  given  by  (4.3 l)i  and  (4.33),  respectively.  Note  that  the  only  difference  in 
this  case  between  the  lower  bound  and  the  upper  estimate  is  that  the  outermost  minimum  in  (4.32) 
is  replaced  by  a  maximum  in  (4.51).  Thus,  from  a  computational  point  of  view,  the  lower  bound 
and  upper  estimates  should  be  computed  together  since  they  involve  the  same  "branches"  which 
must  be  minimized  in  one  case,  and  maximized  in  the  other.  A  special  case  of  the  upper  estimate 
(4.51)  has  been  discussed  previously  by  PONTE  CASTANEDA  (1991b)  for  two-phase  composites 
with  a  linear  and  a  nonlinear  phase.  A  special  case  is  the  rigidly  reinforced  material,  which  has 
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(4.52)  = 

where  phase  2  has  been  selected  as  the  reinforcement. 

Finally,  an  upper  estimate  for  the  incompressible  fiber-reinforced  composite  with 
transversely  isotropic  symmetry  is  given  by  the  same  relation  (4.36)  for  the  lower  bound,  together 
with  (4.37),  except  that  the  minimum  over  all  phases  is  replaced  by  a  maximum. 

4.2.3.  HiRhur-order  upper  estimates.  Higher-order  upper  estimates  for  two-phase  nonlinear 
composites  are  obtained  by  making  use  of  the  lower  bounds  of  Milton  for  comparison  shear 
moduli  with  the  same  distribution  of  phases  as  the  nonlinear  composite  into  (4.7).  These  bounds 
for  the  shear  modulus  and  bulk  modulus  admit  the  same  forms  (4.38)  and  (4.39)  as 
the  corresponding  upper  bounds  and  with  the  difference  that  must  be  replaced  by 
=  2/6.  where  is  given  by  expression  (3.29)  in  MILTON  (1982),  and  is  to  be  replaced  by 

.  These  results  depend  also  only  on  the  two  geometric  parameters  and 
in  addition  to  the  phase  volume  fractions  c-”  and 

As  was  the  case  for  the  lower  bound  the  expression  for  the  corresponding  upper 

estimate  is  very  complicated,  except  for  the  case  where  the  two  phases  arc  incompressible. 
Then,  the  expression  for  reduces  to 


(4.53) 


which  depends  also  only  on  and  we  obtain  the  upper  estimate 
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^(l  +  c^^^cof  +  ^c^^'*rf'^(o^o, 

which  involves  only  one  optimization.  The  corresponding  nonlinear  Hashin-Shtrikman  upper 
estimate  is  obtained  by  choosing  either  =  0,  or  =1,  whichever  yields  the  higher  value. 

4.3.  Attainability  of  the  bounds 

In  this  section,  we  address  briefly  the  question  of  optimality  of  the  nonlinear  Voigt/Reuss 
and  Hashin-Shtrikman  bounds.  By  this  we  mean  whether  specific  microstructures  exist,  within  a 
given  class  of  composites,  which  have  effective  energy  functions  U  attaining  the  bounds  (for  the 
given  class  of  composites).  Further,  in  the  case  of  anisotropic  composites,  microstructures 
saturating  the  bounds  must  exist  for  every  possible  loading  configuration,  in  order  for  the  bound 
to  be  optimal.  If  the  bound  is  saturated  only  for  special  loading  configurations,  the  bound  is  said  to 
be  sharp  (but  not  optimal).  We  remark  that,  even  in  the  context  of  linear-elastic  composites,  there 
are  still  many  open  questions  in  connection  with  the  optimality  of  bounds,  and  most  of  the  results 
that  are  known  to  date  are  for  two-phase  systems.  For  this  reason,  we  will  limit  our  attention  in 
this  section  to  two-phase  nonlinear  composites.  We  further  note  that  in  the  context  of  linear 
elasticity,  the  knowledge  of  optimal  bounds  for  the  effective  energy  functions  of  a  given  class  of 
composites  does  not  suffice  in  general  to  characterize  the  set  of  all  possible  elasticity  tensors  within 
the  class  of  composites,  also  known  as  the  G-closure.  However,  in  some  special  cases,  with 
sufficient  symmetry,  optimal  bounds  for  the  effective  energy  function  do  suffice  to  characterize  the 
G-closure. 

By  construction,  the  nonlinear  bounds  can  only  be  optimal  (sharp),  if  the  corresponding 
linear  bounds  are  optimal  (sharp).  The  microstructures  that  are  commonly  used  in  the  context  of 
linear-elastic  composites  to  prove  optimality  are  sequentially  laminated  materials  (see  KOHN  and 
Milton,  1988).  These  are  hierarchical  microstructures  obtained  by  layering  the  result  of  a 
previous  lamination  procedure  with  one  of  the  original  (or  a  new)  phases  in  a  different  direction. 


(4.41) 


(2)a(2) 
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Thus,  a  rank  I  laminate  (Figure  la)  with  orientation  n,  is  produced  by  layering  the  two  phases  in 
prescribed  proportions,  with  a  certain  layer  thickness  <5/  «  1.  A  rank  II  laminate  (Figure  lb)  is 
obtained  by  layering  the  rank  I  laminate  with  one  of  the  initial  phases  (or  with  a  third  phase)  in  a 
different  direction  n„,  in  given  proportions,  with  layer  thickness  5^  «  1,  such  that  S/ZS,,  «  1. 
Higher  rank  laminates  may  be  obtained  by  iterating  this  procedure.  The  main  advantage  of  this 
class  of  microstructures  is  that  their  effective  properties  can  be  computed  exactly,  because  the  fields 
are  piecewise  constant  within  the  composite. 


Figure  1.  (a)  A  rank-I  laminate  and  (b)  a  rank-II  laminate 


For  the  same  reason,  sequentially  laminated  microstructures  would  also  appear  to  be  useful  in 
the  context  of  nonlinear  composites.  Further,  the  new  variational  principle  seems  to  be  ideally 
suited  for  these  microstructures  also  because  the  fields  are  piecewise  constant,  and  therefore  a 
piecewise  constant  choice  for  the  comparison  moduli  is  bound  to  yield  exact  results.  Thus,  we  can 
see  that  the  effective  energy  0  of  a  two-phase,  nonlinear,  rank-I  laminate  (Figure  la)  can  be 
computed  exactly  via  the  new  variational  principle  in  terms  of  a  two-phase,  linear,  rank  I  laminate. 
However,  for  a  fwo-phase,  nonlinear,  rank  II  laminate  (Figure  lb),  the  effective  energy  of  a  three- 
phase,  linear,  rank  n  laminate  is  required  in  order  to  compute  the  effective  energy  of  the  nonlinear 
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laminate  exactly  in  terms  of  the  new  variational  principle.  (This  is  because  the  fields  are  different 
constants  within  different  layers  in  the  laminate,  even  if  they  are  made  of  the  same  material).  This 
observation  has  the  implication  that  nonlinear  bounds  that  are  obtained  from  optimal  (sharp)  linear 
bounds  involving  sequentially  laminated  microstructures  of  rank  higher  than  one  are  probably  not 
optimal.  This  is  because  the  microstructures  that  attain  the  corresponding  linear  bounds,  from 
which  the  nonlinear  bounds  are  derived,  cannot  be  optimal  in  general  for  the  nonlinear  materials. 
Only  in  the  (unusual)  situation  where  the  optimal  (sharp)  microgeometry  is  a  rank  I  laminate  can 
we  be  sure  of  the  corresponding  optimality  (sharpness)  of  the  nonlinear  bound.  Additionally,  in  the 
special  case  where  the  "matrix"  phase  is  linear,  the  nonlinear  lower  bound  (on  [/ )  is  also  optimal, 
even  when  the  optimal  microstructures  are  laminates  of  rank  higher  than  one  (PONTE  CaSTANEDA, 
1991b).  In  the  next  subsections,  we  will  address  the  optimality  of  the  nonlinear  Voigt/Reuss 
bounds,  and  Hashin-Shtrikman  bounds  by  means  of  the  above-mentioned  criteria. 

4.3.1.  Voigt/Reuss  bounds.  For  linear-elastic,  two-phase,  anisotropic  composites  with 
prescribed  volume  fractions,  the  Voigt/Reuss  bounds  on  the  effective  energy  function  0  are 
known  to  be  non-optimal.  AVELLANEDA  (1987)  has  determined  optimal  bounds  for  the  effective 
energy  function  V  for  this  class  of  composites.  KOHN  and  LiPTON  (1988)  have  found  optimal 
bounds  for  the  subclass  of  incompressible  two-phase  composites,  and  determined  additionally  that 
the  Voigt  (but  not  the  Reuss)  bound  is  actually  optimal  in  this  special  case.  (Although  the  new 
variational  principle  could  in  principle  be  used  in  conjunction  with  these  optimal  bounds  to  obtain 
corresponding  nonlinear  bounds,  these  optimal  linear  bounds  have  rather  complicated  form,  and 
that  is  why  the  simpler  Voigt/Reuss  bounds  were  used  instead).  Thus,  the  nonlinear  Voigt/Reuss 
bounds  are  not  optimal  in  general.  However,  the  linear  Voigt  bound  for  the  two-phase, 
incompressible  composite  does  involve  rank  I  laminate  optimal  microstructure,  and  by  the  criteria 
established  in  the  previous  paragraph,  the  nonlinear  Voigt  bound  is  optimal  for  incompressible 
two-phase  composites  (but  not  more  generall> ).  Although  the  nonlinear  Reuss  bound  is  not 
optimal,  it  is  sharp  under  special  loading  configurations  (see  KOHN  and  LiPTON,  1988  for  the 
corresponding  result  for  linear  elastic  composites).  We  add  that  simple  expressions  for  the  effective 
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properties  of  nonlinear,  rank  I  laminates  have  been  determined  by  DE  BOTTON  and  PONTE 
CaSTASeda  (1991).  Their  results  are  in  agreement  with  the  conclusions  of  this  subsection. 

4.3.1.  Hashin-Shtrikman  lower  bound.  For  linear-elastic,  two-phase,  isotropic  composites 
with  prescribed  volume  fractions,  the  Hashin-Shtrikman  bounds  on  the  effective  energy  function 
0  are  optimal.  This  was  demonstrated  by  Francfort  and  MURAT  (1986).  However,  the 
corresponding  optimal  microstructures  involve  high  rank  laminates,  and  therefore  these 
microstructures  do  not  serve  to  attain  the  nonlinear  Hashin-Shtrikman  lower  bound  on  U  (as 
explained  above).  An  exception  is  provided  by  the  two-phase  composite  with  one  linear  phase;  in 
this  case,  the  lower  bound  can  be  shown  to  be  optimal  (PONTE  CaSTAJ^DA,  1991b).  Finally,  we 
note  that  in  general  the  nonlinear  upper  estimate  for  U  is  probably  as  far  below  the  optimal  upper 
bound,  as  the  lower  bound  is  below  the  optimal  lower  bound  (and  hence  cannot  be  an  upper 
bound).  Thus,  if  we  ignore  the  rigorous  interpret  .tion  of  the  Hashin-Shtriman  estimates  as 
bounds,  and  think  of  them  instead  as  estimates  corresponding  to  the  strongest  and  weakest 
composites,  then  in  practice  both  of  them  should  be  held  on  an  equal  footing.  Thus,  it  appears  that 
it  is  only  on  a  rigorous  sense  that  the  lower  bound  is  "superior"  to  the  upper  estimate. 

5.  Concluding  Remarks 

In  this  paper,  we  have  introduced  dual  versions,  (3.21)  and  (3.34),  of  a  new  variational 
principle  that  allows  the  estimation  of  the  effective  energy  functions  of  nonlinear  composites  with 
isotropic  phases,  as  defined  by  the  classical  minimum  energy  and  minimum  complementary-energy 
principles,  (2.4)  and  (2.11),  respectively.  The  estimation  of  the  effective  properties  is  in  terms  of 
the  effective  energy  functions  of  suitably  optimized  linear  comparison  composites.  These  new 
variational  principles  generalize  a  procedure  introduced  by  PONTE  CASTANEDA  (1991a)  for 
bounding  the  effective  properties  of  nonlinear  composites. 

The  variational  principles  have  been  used  to  derive  bounds  and  estimates  for  the  effective 
properties  of  certain  classes  of  nonlinear  composites.  Thus,  upper  and  lower  bounds  of  the  Voigt 
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and  Reuss  type  have  been  derived  for  the  effective  complementary-energy  function  U  of  generally 
anisotropic  composites,  with  multiple  phases  in  prescribed  volume  fractions.  Although  these 
bounds  are  classical,  new  forms  (4.15)  and  (4.17)  have  been  proposed  for  the  Voigt  bound.  The 
nonlinear  Voigt/Reuss  bounds  were  found  to  be  in  general  not  optimal.  Additionally,  nonlinear 
lower  bounds  of  the  Hashin-Shtrikman  variety  have  been  given  for  the  effective  complementary- 
energy  function  0  of  the  classes  of  isotropic,  and  transversely  isotropic,  fiber-reinforced 
composites  with  prescribed  volume  fractions.  Although  bounds  of  this  type  are  not  new,  and  can 
be  alternatively  obtained  via  the  Talbot-Willis  variational  procedure,  the  form  of  the  bounds  (4.30) 
to  (4.37)  given  here  is  new  and  is  the  simplest  available  thus  far  for  general  multiple-phase 
composites.  We  note  in  this  context  that  most  of  the  explicit  results  to  date  have  dealt  with  porous 
and  rigidly  reinforced  materials.  The  Hashin-Shtrikman  lower  bounds  can  be  shown  to  be  optimal 
(Ponte  Castaneda,  1991b)  for  two-phase  materials  if  one  phase  is  linear,  but  are  probably  not 
optimal  more  generally.  We  have  also  proposed  that  estimates  obtained  by  applying  the  Hashin- 
Shtrikman  upper  bounds  for  the  linear  comparison  composite  into  the  approximate  version  of  the 
new  variational  principle  should  be  interpreted  as  "upper  estimates"  for  the  n  inear  composite, 
corresponding  to  microstructures  exhibiting  the  least  stiff  phase  as  the  matrix  phase.  Finally,  the 
first  nonlinear,  third-order  bounds  of  the  Reran  type  have  also  been  determined  for  the  class  of 
isotropic,  two-phase  composites  with  prescribed  volume  fractions  and  given  Milton  geometric 
parameters.  When  both  phases  are  incompressible,  remarkably  simple  forms  result,  (4.41)  for  the 
lower  bound  and  (4.54)  for  the  upper  estimate,  and  they  reduces  to  the  corresponding  nonlinear 
Hashin-Shtrikman  lower  bound  and  upper  estimate  for  an  appropriate  choice  of  the  Milton 
parameters.  Clearly,  many  possibilities  exist  for  further  investigation.  Central,  and  probably  most 
difficult,  among  them  is  the  determination  of  rigorous  upper  bounds  of  the  Hashin-Shtrikman  and 
higher-order  types  (for  the  effective  complementary-energy  function  of  nonlinear  composites). 
Also,  extensions  of  the  new  variational  principles  are  needed  to  deal  with  nonlinear  composites 
with  anisotropic  phases,  such  as  are  present  in  ductile  polycrystalline  media. 
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Appendix  1. 

The  nonnegative  functions  ^and /introduced  in  §  3  are  defined  on  QxR*  xR*.  On  the  other 
hand,  the  usual  definition  of  the  polar  function  refers  to  functions  defined  on  multiple  copies  of  R 
(the  extended  set  of  reals).  In  this  appendix,  we  show  how  to  define  polars  for  functions  of  the 
type  of  ^  and  f,  defined  on  multiple  copies  of  R* .  For  simplicity,  we  cany  out  the  development  for 
a  one-dimensional  function  (and  drop  the  dependence  on  x);  similar  developments  hold  for 
multiple-dimensional  functions.  Also,  we  consider  only  the  case  of  the  convex  polar  function; 
analogous  results  hold  for  the  concave  polar  function.  Thus,  we  consider  a  function  t:R*  — >  R* 
satisfying  the  condition  t(0)  =  0,  and  define  its  extension  to  R  by,  say 

(1.1)  7(x)=t(|x1). 

where  f;R  R*  is  also  seen  to  be  a  nonnegative  function  satisfying  f(0)  =  0. 

The  convex  polar  function  of  t  is  then  defined  in  the  usual  way  by 

(1.2)  f‘(y)  =  sup{xy-f(x)}. 

Now,  for  y  ^  0,  we  have  that 

f*(y)  =  su^{xy-T(x)} 

(1.3)  =sup{x,.-r(Ar)}_ 

where  we  have  also  used  the  fact  that  t  is  nonnegative  and  t(0)  =  0.  Note  further  that  f*  is  also 
nonnegative,  and  such  that  f*(0)  =  0.  Therefore,  it  makes  sense  to  define  t*:R*  — >  R*  via 

(1.4)  t’(>')  =  sup{xy -  T(x)}. 

Finally,  it  follows  from  the  corresponding  duality  relation  for  f  that 

(1.5)  T(x)>su|»{;0’-'f*(>’)}. 


with  equality  if  t  is  convex. 
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Appendix  II 


The  main  result  proved  in  this  appendix  is  that  concavity  off  implies  convexity  of  g.  The 
demonstration  of  this  result  has  two  parts.  We  begin  by  recalling  some  result  from  §  3.  Thus, 
relation  (3.3)  can  be  written  (dropping  the  dependence  on  x) 


(n.i) 


and  rewritten 


(n.2) 


=  ,sup^o{f 

«(v.,v„)  =  ^sup^„||V^  + 


in  terms  of/ and  g.  Similarly,  relation  (3.6)  can  be  written 


(n.3) 


which  in  turn  yielc 


(n.4) 


=  ^sup^^||a,£,  +  3tT„e„  - 


The  first  part  of  the  demonstration  follows  from 


(n.5) 


=  3VIiX  -«(v.,vj|j 


where  we  have  used  (II.4),  the  fact  that  -sup{/i}=inf{-/i),  and  interchanged  the  resulting  infimum 
over  the  v  variables  with  the  infimum  over  the  u  variables.  But  the  innermost  infimum  is  trivial  to 


compute,  and  there  results 


(11.6) 


p.  4p„  '] 
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which  leads  to  the  important  identity 


(n.7) 


"_1 _ 


The  second  and  final  part  of  the  demonstration  makes  use  of  (II.2)  and  of  the  fact  that  /is 
concave  (so  that  /=  /„)  to  obtain 


Vw  +  3V«r  -  oKp.  +  “-P-  -  Mp-P.)}] 

Viw + -  ".p.  -  “.p.} + /.(p..p-)}. 


where  we  have  'nterchanged  the  suprema.  The  innermost  supremum  is  the  negative  of  the 
innermost  infimum  in  (II.5),  and  therefore  making  use  of  (11.7),  we  deduce  that 


(n.7) 


=  sup  + 
p..p.2o[9p,  4p^ 

=  5*'K.v„). 


1  _9_ 


This  final  result  proves  that  g  is  convex,  and  hence  completes  the  proof. 
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Appendix  III 


In  this  appendix,  we  give  a  simple  proof  of  the  identity  (4.12)  used  throughout  the  body  of  the 
paper.  We  begin  by  introducing  a  Lagrange  multiplier  to  account  for  the  zero-average  constraint  on 
the  optimization  variables,  CD^'^.  Thus, 


and,  by  means  of  the  Saddle  Point  theorem,  we  have 


<usO  V 


c‘'W'’ 


The  infimum  over  the  O)^'^  is  satisfied  by 


(III.3) 


which  finally  leads  to 


2 


s=0  t  '■=1 


i!?; 


(III.4) 


Q  (O  ( 

=  sup /I  -  —  ^  c 


Reference  [12] 


In  Anisotropy  and  Localization  of  Plas^c  Deformadon,  edited  by  J.-P.  Boehler 
and  A.  S.  Khan,  pp.  15-18.  London;  Elsevier,  1991. 
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ABSTRACT 

A  recently  developed  variational  principle  for  estimating  the  effective  properties  of  nonlinear 
composites  in  terms  of  the  corresponding  properties  of  linear  composites  with  the  same 
microstructural  distributions  of  phases  is  applied  to  two  model  anisotropic  composite  materials. 
The  model  materials  considered  are  laminat^  materials  and  fibre-reinforced  materials,  and  they 
correspond  in  the  dilute  limit  of  the  inclusion  phase  to  materials  reinforced  by  aligned  platelets 
and  fibers,  respectively.  For  simplicity,  the  power  exponent  of  one  of  the  phases  will  be  taken 
to  be  unity  correspondng  to  linear  behavior,  while  Ae  other  phase  will  assumed  to  satisfy  a 
pure-power  creeping  law.  Both  phases  will  be  assumed  to  be  isotropic  and  incompressible. 


INTRODUCTION 

A  new  procedure  for  estimating  the  effective  properties  of  composite  materials  with  phases 
exhibiting  nonlinear  constitutive  behavior  has  been  proposed  recently  by  Ponte  Castaneda  [6]. 
The  straightforward  and  versatile  procedure  expresses  fte  effective  properties  of  the  nonlinear 
composite  in  terms  of  the  effective  properties  of  a  family  of  linear  composites  with  the  same 
distribution  of  phases  as  the  nonlinear  composite.  Thus,  bounds  and  estimates  for  the  effective 
properties  of  linear  composites  can  be  translated  directly  into  bounds  and  estimates  for  the 
corresponding  nonlinear  composite.  Appropriate  references  for  the  linear  theory  of  composites 
are  provided  by  the  works  of  Christensen  [1].  The  new  procedure  was  applied  in  the  above 
reference  [6]  to  composite  materials  containing  a  nonlinear  isotropic  matrix  either  weakened  or 
reinforced  by  isotropic  distributions  of  voids  or  rigid  particles,  respectively.  The  case  of  a 
ductile  matrix  reinforced  by  incompressible  elastic  particles  was  studied  in  [7].  The  results  of 
these  studies  were  given  in  the  form  of  estimates  and  rigorous  bounds  for  the  effective 
properties  of  such  materials.  The  Hashin-Shtrikman  bounds  obtained  via  the  new  method 
directly  from  the  Hashin-Shtrikman  [3]  bouni* .  for  the  linear  comparison  material  were  found 
to  be  an  improvement  over  the  corresponding  bounds  obtained  in  [8]  for  the  same  class  of 
nonlinear  materials  using  the  nonlinear  extension  of  the  Hashin-Shtrikman  variational  principle 
[3]  proposed  by  Talbot  and  Willis  [9].  Recendy,  however,  Willis  [11]  has  shown  that  the 
bounds  obtained  via  the  new  method  can  also  be  obtained  by  the  TallMt-Willis  method  with  a 
better  choice  of  their  comparison  material.  More  generally,  the  new  procedure  can  make  use  of 
other  bounds  and  estimates  for  the  linear  comparison  material  to  yield  corresponding  bounds 
and  estimates  for  nonlinear  materials.  In  fact,  the  new  procedure  can  be  shown  [2]  to  yield 
exact  results  for  nonlinear  composites  with  speci microstructures. 

In  this  paper,  we  apply  the  new  procedur  composite  materials  containing  a  creeping 
phase  of  the  power-law  type  and  a  linearly  creeping  phase.  Although  the  phases  will  be 
assumed  to  be  isotropic,  the  composite  itself  will  be  assumed  to  be  anisotropic  by  identifying 


preferred  orientations  in  the  microstructural  distribution  of  the  phases.  Two  examples  will  be 
considered,  both  of  which  have  transversely  isotropic  overall  symmetries:  laminat^  materials 
and  fibre-reinforced  materials.  The  first  class  of  materials  can  be  given  the  interpretation  of 
materials  reinforced  by  aligned  platelets  for  dilute  concentrations  of  the  linear  phase;  whereas, 
the  second  class  of  materials  can  be  given  the  interpretation  of  materials  reinforced  by  short 
fibers  for  dilute  concentrations  of  the  linear  phase.  For  simplicity,  we  assume  that  the  phases 
are  incompressible  and  perfectly  bonded  to  each  other.  In  the  next  three  sections,  we  will  give, 
respectively,  a  brief  definition  of  effective  properties,  an  introduction  to  the  new  variational 
principle  of  [6],  and  application  of  the  principle  to  the  two  classes  of  materials  discussed 
above.  More  detailed  studies  including  specific  results  for  the  nonlinear  laminates  and  fiber- 
reinforced  materials  can  be  found  in  [2]  and  [5],  respectively. 


EFFECTIVE  PROPERTIES 

Consider  a  two-phase  composite  occupying  a  region  of  unit  volume  A  such  that  the  local 
stress  potential  U(a,x)  is  expressed  in  terms  of  the  homogeneous  phase  potentials  via 

U(a,x)  =  ^x''\x)U^'\a),  (1) 

r=l 

where  is  the  characteristic  function  of  phase  r.  The  phases  are  assumed  to  be 

incompressible  and  isotropic,  so  that  the  potentials  U^'\a)  can  be  assumed  to  depend  only  on 

the  effective  stress  o;  =  5 ,  where  S  is  the  deviator  of  a .  Thus,  we  write 

where  the  are  scalar-valued  functions.  Then,  the  local  constitutive 
relation  for  the  creeping  material  is  given  by 

D  =  -^(a,x),  (2) 

d<5 

where  D  is  the  rate-of-deformation  (strain-rate)  tensor. 

To  define  the  effective  properties  of  the  heterogeneous  material  we  introduce,  following 
Hill  [4],  a  uniform  constraint  boundary  condition 

on  =  an,  X€(3Q,  (3) 

where  dQ.  denotes  the  boundary  of  the  composite,  n  is  its  unit  outward  normal,  and  o  is  a 
given  constant  symmetric  tensor.  It  follows  that  the  average  stress  is  precisely  a,  and  we 

drfine  the  average  strain-rate  D  in  a  similar  manner. 

Then,the  effective  behavior  of  the  composite,  or  the  relation  between  the  average  stress  and 
the  average  strain-rate  follows  from  the  principle  of  minimum  complementary  energy,  which 
can  be  stated  in  the  form 


{7(a)=  min  \u{cs,x)dV,  (4) 

a€5(<j)JQ 

where  5(c)  =  {alcr^  ^  =0  in  A  and  a^rij  =  a^rij  on  is  the  set  of  statically  admissible 
stresses,  and  where  we  have  assumed  convexity  of  the  nonlinear  potential  C/(a,x).  Thus, 
assuming  that  U(o)  is  differentiable,  we  have  that 


D  =  -p(a,x). 
da 


(5) 


The  task  will  be  to  determine  bounds  and  estimates  for  (/(a),  which  is  known  to  be  convex. 


NONLINEAR  VARIATIONAL  PRINCIPLES 


A  new  variational  principle  for  determining  bounds  and  estimates  for  the  effective  properties  of 
nonlinear  composites  in  terms  of  the  effective  properties  of  linear  composites  was  proposed  by 
Ponte  Castaneda  [6].  In  this  section,  we  specialize  this  result  for  the  case  where  bo  A  phases 
are  incompressible,  and  phase  2  is  linear  so  that 


The  new  variational  principle  is  based  on  a  representation  of  the  potential  of  the  nonlinear 
material  in  terms  of  the  potentials  of  a  family  of  linear  comparison  materials.  Thus,  for  a 


homogeneous  nonlinear  material  with  "stronger  than  quadratic"  growth  in  its  potential,  U(g), 
and  certain  additional  convexity  hypothesis,  we  have  that 


where 


C/(a)  =  max{f/,(a)-V(//)}  , 

ft>0 


V(M)  =  tnax[U,(a)-Uia)] 

O 


(6) 

(7) 


and  where  U^(a)  is  the  potential  of  a  linear  comparison  material  with  shear  modulus  p . 

The  new  variational  principle  is  obtained  by  making  use  of  relation  (6)  applied  to  the 
nonlinear  phase  1  in  the  complementary  energy  principle  (4)  to  obtain  the  following  relation  for 
the  nonlinear  composite 


f7(o)  =  nipfe(CT)-  /V<»(At»>)dv|, 

where 


(8) 


l7,(a)=  nun  f  t/,(a,x)dV,  (9) 

is  the  effective  potential  of  a  linear  comparison  material  with  local  potential  f/„(a,x)  and  shear 
moduli  and  in  phases  1  and  2,  respectively.  Further,  we  note  that,  in  general,  the 
comparison  moduli  are  functions  of  x. 

The  variational  principle  described  by  (8)  roughly  corresponds  to  solving  a  linear  problem 
for  a  heterogeneous  material  with  arbitrary  moduli  variation  within  the  nonlinear  phase,  and 
then  optimizing  with  respect  to  the  variations  in  moduli  within  the  nonlinear  phase.  Thus,  the 
nonlinear  material  can  be  thought  of  as  a  "linear"  material  with  variable  moduli  that  are 
determined  by  prescription  (8)  in  such  a  way  that  its  properties  agree  at  each  x  with  those  of  the 
nonlinear  material,  lliis  suggests  that  if  the  Helds  happen  to  be  constant  over  the  nonlinear 

phase,  then  the  variable  moduli  can  be  replaced  by  constant  moduli  p^'\  More 

generally,  however,  we  have  the  following  lower  bound  for  U{a),  namely, 

(7_(a)  =  max{17,(a)-c“V"'(/i''')},  aO) 

where  is  the  volume  fraction  of  phase  1. 


APPLICATION  TO  LAMINATES 


In  this  section,  we  summarize  from  reference  [2]  the  results  of  applying  the  new  variational 
principle  to  a  laminated  material.  For  this  geometry,  it  is  well  known  that  the  stress  and  strain 

fields  take  on  different  constant  values  within  each  phase.  This  suggests  that  taking  to  be 
constant  within  the  nonlinear  phase  will  lead  to  an  exact  result  as  given  by  (8)  or  (10)  (in  this 

case  U  and  U_  are  identical).  The  key  ingredient  in  this  development  is  the  exact  result  for  the 
effective  potential  of  the  linear  laminate.  This  result  can  be  expressed  in  the  form 

(11) 

where  T„,  T^,  o;  and  are  the  transversely  isotropic  invariants  of  the  stress  tensor 
corresponding  respectively  to  the  out-of-plane  tensile,  in-plane  hydrostatic,  out-of-plane  shear 

and  in-plane  shear  stresses.  Also,  fi  and  fi  stand,  respectively,  for  the  Voigt  and  Reuss 
estimates  of  the  shear  modulus. 

Application  of  this  result  into  the  new  variational  principle  can  then  be  shown  to  lead  to  the 
following  simple  expression  for  the  effective  energy  of  the  laminate 

UCc)  =  (12) 


where  s”'  and  are  functions  of  co  given  by  the  relations  =  -y (l  + 
and  =  -^(l-c‘”<u)^(<T^  “  'r^)+  • 

Analogous  forms  can  be  derived  for  the  fibre-reinforced  material,  except  that  for  this 
geometry,  the  effective  energy  cannot  be  given  explicitly  by  (8),  but  only  characterized  in  the 
form  of  bounds  by  means  of  (10).  Detailed  results  of  these  calculations  are  given  in  [5],  where 
they  are  also  compared  to  the  results  of  Talbot  and  Willis  [10]  for  the  same  microstructure 
using  the  Talbot-Willis  variational  principle.  Finally,  both  types  of  results  can  be  averaged  over 
all  possible  orientations  to  yield  results  depicting  approximately  the  effect  of  inclusion  shape  on 
the  effective  properties  of  isotropic  composites. 
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Abstract 

A  laminated  material  is  one  of  the  few  composite  systems  for  which  the  effective  constitutive 
behavior  can  be  determined  exactly.  This  is  well-known  for  laminated  composites  with  linear-elastic 
phases  in  prescribed  volume  fractions,  for  which  explicit  expressions  have  been  available  for  their 
effective  moduli  for  at  least  thirty  years.  However,  it  appears  that  corresponding  expressions  for  the 
effective  energy  functions  of  laminated  composites  with  phases  exhibiting  nonlinear  constitutive 
behavior  are  currently  unavailable.  While,  in  principle,  the  calculation  of  the  effective  energy 
function  of  a  laminated  composite  is  straightforward;  in  practice,  the  standard  procedure  involves  the 
solution  of  a  fairly  complicated  set  of  nonlinear  algebraic  equations.  In  this  paper,  we  make  use  of  a 
new  variational  procedure  recently  developed  by  one  of  the  authors  to  obtain  fairly  explicit 
expressions  for  the  effective  energy  functions  of  laminated  composites  with  ductile  phases  in 
prescribed  volume  fractions.  Explicit  results  are  then  computed  for  ductile/brittle  systems,  such  as  an 
aluminum/alumina  laminate,  and  for  a  laminated  composite  made  up  of  two  rigid/perfectly  plastic 
phases  with  different  yield  stresses.  The  results — ^which  are  representative  of  other  anisotropic 
geometries,  such  as  fiber-reinforced  solids — exhibit  a  strong  anisotropic  coupling  between  different 
loading  modes  that  is  enhanced  by  material  nonlinearity. 
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1.  Introduction 


This  work  is  concerned  with  the  determination  of  the  effective  behavior  of  laminated  composite 
materials  with  plastically  deforming  phases  in  prescribed  volume  fractions.  We  will  deal  with  the 
exact  effective  behavior  of  such  materials,  and  therefore,  we  will  exclude  from  our  consideration  the 
so-called  approximate  theories  of  laminated  plates  (Christensen  1979).  In  the  context  of  linear 
elasticity,  laminated  composites  have  been  studied  extensively  in  the  literature.  Postma  (1955)  and 
White  and  Angona  (1955)  concerned  themselves  with  the  study  of  two-phase,  periodic  imninates  in 
connection  with  wave  propagation  in  stratified  media.  Backus  (1962)  extended  their  results  to  multi¬ 
phase,  non-periodic  composites,  again  in  the  wave  propagation  context.  The  extension  to  anisotropic 
layers  was  considered  by  Walpole  (1969)  for  aligned,  transversely  isotropic  phases,  and  by  Chou  et 
al.  (1972)  and  Pagano  (1974)  for  more  general  phase  anisotropies.  Recently,  Norris  (1991)  has 
developed  alternative  expressions  for  the  effective  moduli  tensor  of  laminated  composites  with 
generally  anisotropic  phases,  exploiting  the  interior  and  exterior  projection  tensors  of  Hill  (1972, 
1983).  Other  related  works  include  iterative  formulae  developed  by  Francfort  and  Murat  (1986)  in 
the  context  of  linear  elasticity  allowing  simple  expressions  for  the  effective  moduli  of  multi-sandwich 
structures  O^UTunates  embedded  within  laminates  of  different  orientations).  These  microstructures  are 
of  theoretical  value  in  the  demonstration  of  the  optimality  of  bounds  for  the  effective  properties  of 
composite  materials  with  more  general  microstructures  (see  Kohn  1987,  Lipton  1991b).  In  spite  of 
the  large  level  of  activity  for  linear  laminates,  briefly  sununarized  above,  the  theory  of  nonlinear 
laminated  composites  does  not  seem  to  have  been  developed  nearly  as  much.  To  the  knowledge  of 
the  authors,  the  only  work  thus  far  in  this  direction  is  a  generalization  of  the  Franefort-Murat  formula 
for  two-phase  laminated  materials  by  Kohn  (1990)  and  Milton  (1990).  Nonlinear  results  do  exist, 
however,  within  the  context  of  laminated  plate  theories. 

The  justification  of  the  study  of  nonlinear  laminated  composites  may  be  partially  understood  in 
terms  of  the  following  considerations.  First,  it  is  a  configuration  of  practical  importance:  for 
instance,  the  use  of  linear  laminated  theories  in  geophysical  applications  is  well-known,  but  it  is  also 
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known  that  the  properties  of  the  materials  composing  the  surface  of  the  Earth  may  also  exhibit 
nonlinear  constitutive  behavior,  particularly,  deep  within  the  surface,  where  the  materials  are  subject 
to  large  compressive  stresses.  Second,  the  laminated  microstructure  corresponds  to  the  simplest 
possible  type  of  anisotropic  composite  with  nonlinear  phases,  in  the  sense  that  exact  results  may  be 
obtained  for  its  effective  properties,  as  will  be  shown  herein.  In  this  connection,  the  present  work 
should  be  considered  in  the  light  of  a  research  program  attempting  to  characterize  the  effective 
properties  of  nonlinear  composites  in  general.  Thus  far,  methods  have  been  developed  for 
understanding  the  effective  behavior  of  nonlinear  composites  by  Talbot  and  Willis  (1985),  and  by 
Ponte  Castaneda  (1991a,  c).  These  methods,  although  different  in  essence,  can  be  shown  to  yield 
exactly  the  same  results  in  some  cases,  but  the  new  method  is  more  general  than  the  fu^t  in  that  it  can 
be  used  to  obtain  estimates  other  than  Hashin-Shtrikman  bounds  and  self-consistent  estimates  (see 
Willis  1991,  and  Ponte  Castaneda  1991a,  c).  These  methods  have  been  applied  to  obtain  bounds  and 
estimates  for  the  effective  properties  of  nonlinear  composites  with  isotropic  overall  synunetries  by 
Ponte  Castaneda  and  Willis  (1988),  and  Willis  (1989)  making  use  of  the  Talbot-Willis  method,  and 
by  Ponte  Castaneda  (1991a,  b,  c)  making  use  of  the  new  method.  Results  for  fiber-reinforced 
composites  have  also  been  developed  very  recently  by  Talbot  and  Willis  (1991),  and  by  Ponte 
Castaneda  (1991c),  and  Ponte  Castaneda  and  de  Botton  (1991).  Third,  we  will  find  that  the 
laminated  microstructure  is  illustrative  of  the  significant  coupling  that  may  arise  in  nonlinear, 
anisotropic  materials  between  different  loading  modes.  This  nonlinear  coupling  is  also  observed  in 
other  more  complex  microstructures,  such  as  fiber-reinforced  composites. 

The  rest  of  the  paper  is  structured  as  follows.  In  §  2,  the  definition  of  effective  properties  is 
reviewed,  and  their  variational  characterization  is  given  in  terms  of  both  the  classical  variational 
principles,  and  the  new  variational  principles  of  Ponte  Castaneda  (1991a,  c).  In  §  3,  the  general 
nonlinear  laminated  composite  is  considered,  and  general  formulae  for  the  effective  properties  are 
derived  in  §  4  and  5  for  incompressible  and  compressible  laminates,  respectively.  Finally,  in  §  6, 
more  specific  results  are  given  for  two-phase  laminates.  In  particular,  the  cases  of  ductile  materials 
reinforced  by  linear-elastic  layers,  and  of  laminates  with  two  perfectly  plastic  phases  are  considered. 
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2.  Effective  properties  and  their  variational  characterization 


In  this  section,  we  are  interested  in  the  characterization  of  the  effective,  or  overall,  constitutive 
behavior  of  composites  materials  with  plastically  deforming  phases.  For  our  purposes,  a  composite 
is  a  heterogeneous  material  with  two  distinct  length  scales:  one  macroscopic,  L,  describing  the  gross 
size  of  the  specimen  and  the  scale  of  variation  of  the  appUed  loading  conditions,  and  a  microscopic 
scale,  /,  characterizing  the  size  of  the  typical  inhomogeneity,  such  that  / «  L.  More  precise 
definitions  can  be  found  in  the  review  article  by  Kohn  (1987). 

For  simplicity,  the  constitutive  behavior  of  the  phases  will  be  characterized  by  the  deformation 
theory  of  plasticity,  or  equivalently  by  nonlinear  infinitesimal  elasticity.  However,  with  minor 
changes  in  notation,  the  results  of  the  analyses  of  this  paper  will  also  be  relevant  to  the  high- 
temperature  creeping  behavior  of  composite  laminates.  Additionally,  the  results  can  be  used  in  an 
approximate  fashion  to  suggest  yield  functions  for  laminated  composites  in  the  context  of  the 
incremental  theory  of  plasticity,  as  suggested  by  Duva  and  Hutchinson  (1984),  and  other 
investigators. 

In  the  following  description  of  effective  properties,  the  composite  is  assumed  to  occupy  a 
domain  of  unit  volume  with  boundary  dQ.  Then,  the  nonlinear  plastic  behavior  of  the  composite 
is  characterized  by  means  of  a  complementary-energy  density  function,  l/(x,a),  depending  on  the 
position  vector  x  and  the  stress  field  0(x),  in  such  a  way  that  the  strain  field  e(x)  is  given  by 


(2.1) 


e(x) 


dU{\,G) 

do 


Following  Hill  (1963),  we  define  the  effective  constitutive  behavior  of  the  heterogeneous 
solid  in  terms  of  the  analogous  relation 


(2.2) 
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where  E  denotes  the  mean  value  of  the  strain  field  over  Cl,  and  0  refers  to  the  normalized  (recall  that 
Q  has  unit  volume)  complementary-energy  function  of  the  solid  when  subjected  to  the  uniform 
constraint  boundary  condition 

(2.3)  an  =  CTn,  xedCl, 

where  ii  is  the  outward  unit  normal  to  dCl,  and  0  is  a  constant,  synunetric  tensor.  We  recall  that 
under  this  type  of  boundary  condition,  the  mean  value  of  the  stress  over  Cl  is  precisely  o. 

The  effective  complementary-energy  function  of  the  composite,  U,  can  be  obtained  directly  in 
terms  of  the  principle  of  minimum  complementary  energy  by  means  of 

(2.4)  U(a)=  min  fl/(x,a)dv, 

G€5(c)  I 
Q 

where 

(2.5)  5(a)  =  {a|  V  0  =  0  infl,  and  Gn  =  OT  on 

is  the  set  of  statically  admissible  stress  fields.  Note  that  the  first  set  of  conditions  in  (2.5)  are  the 
equilibrium  equations,  and  that  the  minimizing  conditions  (Euler-Lagrange  equations)  of  (2.4)  are  the 
compatibility  equations.  Further,  composite  materials  typically  exhibit  sharp  interfaces  across  which 
the  material  properties  are  discontinuous,  although  the  phases  are  assumed  to  be  perfectly  bonded. 
Therefore,  across  such  interfaces,  the  equilibrium  equations  must  be  reinterpreted  in  terms  of 
continuity  of  the  traction  stresses,  and  correspondingly  the  compatibility  equations  must  be  replaced 
by  continuity  of  the  tangential  components  of  the  strain  tensor . 

We  note  that,  given  relation  (2.2)  in  terms  of  f7 ,  the  problem  of  characterizing  the  effective 
behavior  of  the  composite  reduces  to  that  of  determining  U.  However,  while  in  principle  U  can  be 
computed  from  (2.4);  in  practice,  this  variational  principle  is  not  very  useful  for  two  reasons.  First, 
usually  the  microstructure  of  a  typical  composite  is  not  completely  specified;  and  second,  the 
problem  described  by  (2.4)  is  a  nonlinear  one  in  account  of  the  nonlinear  behavior  of  the  constituent 
phases.  For  the  problem  of  interest  in  this  paper,  the  first  issue  is  not  a  concern  because  the  phase 
volume-fractions  suffice  to  characterize  the  microstructure  of  a  laminated  composite  material. 
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However,  the  second  issue  presents  real  difficulties.  For  this  reason,  we  describe  next  a  new 
variational  principle,  introduced  recently  by  Ponte  Castaneda  (1991a,  c),  which  deals  precisely  with 
the  problem  of  constitutive  nonlinearity.  This  is  accomplished  by  expressing  the  effective  energy 
function  of  the  nonlinear  composite  in  terms  of  a  variational  statement  involving  the  effective  energy 
functions  of  the  class  of  linear  comparison  composites.  Thus,  the  new  variational  principle  allows 
the  extension  of  well-known  results  for  linear  composites  to  corresponding  results  for  nonlinear 
ones.  In  this  paper,  we  will  make  use  of  this  variational  principle,  and  of  well-known  results  for  the 
effective  properties  of  linear-elastic  laminates,  to  determine  the  effective  constitutive  behavior  of 
ductile  laminates.  Before  proceeding  with  this  task,  we  briefly  review  the  new  variational  principle. 

The  new  variational  principle  for  the  effective  energy  of  the  composite  U  is  obtained  by  means 
of  the  Legendre  transformation,  applied  to  a  modified  set  of  variables.  We  will  assume  that  the 
heterogeneous  solid  is  (locally)  isotropic,  such  that 
(2.6)  C/(x,o)=  v'(x;t„<t„), 

where  is  a  nonnegative  function,  that  is  convex  in  its  last  two  arguments,  and  satisfies  the 
condition  that  y/{x;0,0)  =  0  for  all  x ,  and  where 


(2.7) 


=-tra  and 

m  3 


are  the  mean  and  effective  shear  (in  its  plasticity  denotations)  stresses,  respectively,  with 
o'  =  a  -  <T„I  denoting  the  stress  deviator  tensor.  We  note  that  form  (2.6)  is  not  the  most  general 

form  for  the  energy  function  of  a  nonlinear  isotropic  solid  (we  could  also  have  dependence  on  the 
determinant  of  O),  but  this  form  is  still  general  enough  to  cover  the  usual  plasticity  models  of  interest 
here.  Further,  we  will  assume  that  the  growth  in  C/  as  the  magnitude  of  the  stress  tensor  becomes 
large  is  stronger  that  quadratic.  This  is  of  course  consistent  with  the  ductile  behavior  of  the  material. 

The  new  variational  principle  is  obtained  then  in  terms  of  the  following  expression  for  the 
energy-density  function  of  the  heterogeneous  solid,  namely, 

(2.8)  Uix,a)= 
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where  is  the  energy-density  function  of  a  linear-elastic  comparison  solid  with  shear  modulus 
and  bulk  modulus  k^,  such  that 


(2.9)  = 
and  where 

(2.10)  V(x-,n^,Kj  =  m2a{f/^(x,a)-l/(x,o)}. 

(Note  that  the  maximum  in  the  above  function  is  usually  bounded,  because  of  the  stronger  than 
quadratic  assumption  on  U).  These  expressions  are  obtained  by  means  of  the  changes  of  variables, 

V,  =  and  ,  which  lead  to  the  deflnition  of  a  nonnegative  function/,  such  that 

(2.11)  /(x;v„v„)=  v'(x;T„cr„). 

Then,  expression  (2.10)  is  nothing  more  than  the  Legendre  dual  of  f;  in  fact. 


V(x:'t^,x^)  =  /*(x;p„p„),  with  p,  and  p„  =-^,  Here,  /'  is  the  Legendre  transform  of/ 


2/t 


2k ^ 


given  by 

(2- 12)  /’(x; 

and  (2.8)  is  a  statement  of  Legendre  duality /or  convex/ {i.e.,  /“  =  /,  but  written  in  the  terms  of  U 
and  V).  For  details,  we  refer  the  reader  to  Ponte  Clastaneda  (1991c). 

The  new  variational  principle  is  then  obtained  essentially  by  inserting  expression  (2.8)  for  U 
into  the  principle  of  minimum  complementary  energy  (2.4),  and  interchanging  the  order  of  the 
minimum  over  the  set  of  admissible  stresses  with  the  maximum  over  the  comparison  moduli.  The 
result  may  be  expressed  in  the  form  (Ponte  Castaneda  1991c) 


(2.13) 


^(a)  = 


max 


Jv(x;p„(x),j(:„(x))d[x  -, 

Q 
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where 


(2.14)  Lr„(a)=  imn  lf/,(x,a)dlx 

aeS(0)  J 

a 

is  the  effective  energy  of  the  linear  comparison  composite.  We  emphasize  that  expression  (2.13)  is  a 
variational  principle  in  its  own  right  since  it  involves  an  infinite-dimensional  optimization  over  the  set 
of  nonnegative  functions  /io(x)  and  k'^(x).  Thus,  even  if  we  had  an  explicit  expression  for  the 
effective  energy  function  of  the  linear  comparison  composite  (not  an  easy  calculation  in  general 
for  arbitrary  /t„(x)  and  the  above  variational  principle  would  still  be  difficult  to  implement. 

However,  we  will  see  that  for  a  laminated  composite,  the  above  problem  simplifies  dramatically. 
Similarly,  Ponte  Castaneda  (1991a,  b,  c)  have  shown  that  the  above  variational  principle  can  also  be 
utilized  in  an  approximate  fashion  to  compute  bounds  for  the  effective  properties  of  nonlinear 
composites  with  more  general  microstructures.  In  the  same  references,  dual  versions  of  (2.13)  are 
also  given  in  terms  of  the  minimum  potential  energy  of  the  composite;  however,  in  this  paper  we 
prefer  to  use  the  above  formulation  due  to  the  fact  that  it  is  easier  to  express  the  stress/strain  relation 
for  a  ductile  material  in  terms  of  the  complementary  energy-density  function  U  than  in  terms  of  its 
Legendre  counterpart,  the  energy-density  function  W  -  U*. 

3.  Application  of  the  new  variational  principle  to  a  nonlinear  laminated  composite 

In  this  section,  we  specialize  the  general  formulation  of  the  previous  section  to  the  case  of  a 
laminated  composite.  Such  materials  consist  of  n  homogeneous,  isotropic  phases  occupying  non¬ 
intersecting  layered  regions  (r  =  1,  2,...,  n),  with  union  £l  and  with  normal  n.  The 
complementary  energy-density  function  for  the  laminated  material  is  then  expressible  in  the  form 

n 

(3.  U{G,x)  =  ^/'\x  •  n)  , 

r=I 
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where  ‘  **)  (equal  to  1  for  x  in  phase  r,  and  0  otherwise)  is  the  characteristic  function  of  phase 

r,  and  U^'\a)  =  is  the  corresponding  homogeneous,  isotropic  energy-density  function. 

Also  the  volume  fraction  of  each  phase  is  determined  by  the  corresponding  characteristic 
functions  via  the  relation 

(3.2)  ^^  =  ^x^'\x-n)dx. 

a 

We  remark  that  a  laminated  composite  with  perfectly  bonded,  isotropic  phases  possesses 
transversely  isotropic  symmetry  (with  transverse  direction  n).  In  some  sense,  it  represents  the 
simplest  composite  material  with  transverse  isotropy;  other  examples  of  practical  importance  include 
fiber-reinforced  materials  with  isotropic  constituent  phases.  These  will  be  considered  elsewhere. 
Because  of  the  particular  type  of  anisotropy  involved  in  laminated  composites,  we  have  included  in 
Appendix  I  a  brief  summary,  largely  after  Walpole  (1981),  of  the  appropriate  invariants  and  other 
useful  definitions  for  transversely  isotropic  materials. 

The  computation  of  the  effective  energy-density  function  of  a  laminated  composite  is  made  easy 
by  the  following  property  of  laminated  composites.  If  the  thickness  of  the  typical  layer  is  small 
compared  to  the  size  of  the  laminate  (i.e.,  if  the  laminate  is  a  composite  in  the  sense  defined  in  §  2), 
then,  away  from  a  boundary-layer  region  close  to  the  boundary  of  the  composite,  the  fields  are 
constant  within  each  layer  (a  different  constant  within  each  layer).  Therefore,  the  problem  of 
determining  the  effective  energy  function  of  a  laminated  composite  reduces  to  that  of  determining  the 
constant  fields  within  each  phase  of  the  composite  by  imposition  of  the  appropriate  jump  conditions 
(continuity  of  traction  stresses  and  tangential  strains)  across  the  interfaces  between  the  different 
layers,  as  well  as  the  averaging  conditions  stated  in  §  2.  Thus,  the  problem  of  determining  the 
effective  energy  function  of  a  laminated  composite,  unlike  the  corresponding  problem  for  a  general 
composite,  simplifies  to  an  algebraic  one.  Although,  in  principle,  the  problem  can  always  be  solved; 
in  practice,  it  may  be  difficult  to  obtain  explicit  results  because  the  jump  conditions  take  the  form  of 
complicated  sets  of  nonlinear  algebraic  equations.  However,  if  the  composite  is  made  up  of  linear 
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phases  (with  quadratic  energy  functions  in  each  phase),  the  jump  conditions  are  also  linear  and  they 
can  be  solved  in  closed  form,  as  discused  in  §  1. 

The  results  for  the  effective  energy  functions  of  linear  laminated  composites  are  given  in  the 
next  two  sections;  in  the  balance  of  this  section,  we  make  use  of  the  new  variational  principle  (2.13) 
to  determine  an  exact  expression  for  the  effective  energy  function  of  a  nonlinear  laminated  composite 
U ,  written  in  terms  of  the  effective  energy  functions  of  the  class  of  linear-elastic  comparison 
laminates.  This  is  accomplished  by  noting  that  the  minimizing  comparison  moduli  functions  //<,(x) 
and  K’o(x)  in  (2.13)  must  be  constant  within  each  phase.  Therefore,  it  suffices  to  optimize  with 
respect  to  the  set  of  constant  [over  each  phase  r(  r  =  1, ...,  /i)]  comparison  moduli,  and 
The  new  variational  principle  then  leads  to 


(3.3) 


U(a)=  max 


*»1 


where,  from  (2.10), 


and  where  U„  is  the  effective  energy-density  function  of  a  linear-elastic  laminated  material  made  up 
of  n  phases  in  volume  fractions  c^'\  with  shear  and  bulk  moduli,  and  respectively. 

On  the  face  of  it,  expressions  (3.3)  with  (3.4)  for  the  effective  energy  function  of  a  nonlinear 
laminated  composite  do  not  appear  to  offer  much  of  an  analytical  advantage  over  the  standard 
procedure  of  determining  the  stress  fields  within  each  phase  (by  solving  the  appropriate  nonlinear 
jump  conditions)  and  putting  them  directly  in  the  complementary  energy  principle  (2.4).  This  is  due 
to  the  large  number  of  optimizations  involved  in  expressions  (3.3)  and  (3.4)  {i.e.,  a  total  of  4n 
optimizatior  for  an  n-phase  laminate).  However,  we  shall  see  in  the  next  two  sections  that 
application  oi  tlie  particular  form  for  the  effective  energy  function  of  a  linear-elastic  laminate  in  (3.3) 
leads  to  a  simpler  optimization  problem  for  the  effective  energy  function  of  the  nonlinear  laminate. 
Further,  we  observe  that,  from  a  computational  point  of  view,  it  is  generally  easier  to  minimize 
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functions  than  it  is  to  solve  nonlinear  sets  of  equations,  and  therefore,  the  methods  developed  in  this 
paper  are  computationally  superior  to  the  standard  procedure  of  solving  systems  of  nonlinear 
equations.  In  §  4,  we  begin  by  considering  the  simpler  case  of  a  laminated  composite  with 
incompressible,  isotropic  phases,  and  in  §  5,  we  tackle  the  more  complicated  problem  of  a  general 
laminated  composite  with  compressible,  isotropic  phases. 


4.  The  incompressible  laminated  composite 


In  this  section,  we  deal  with  the  special  case  of  laminated  composites  with  incompressible 
phases.  In  this  case,  the  energy-density  functions  of  each  phase  take  the  simpler  form 
U^'\g)  =  Then,  the  relations  (3.3),  (3.4)  expressing  the  effective  energy  function  i7  of 

the  nonlinear  laminate  reduce  to 


(4.1) 


where 


(4.2) 


(7(o)  =  max 


n 


and  where  now  refers  to  the  effective  energy-density  function  of  a  linear-elastic  laminated 
material  composed  of  n  incompressible  phases  in  volume  fractions  with  shear  moduli  . 

The  effective  energy-density  function  of  the  linear  (incompressible)  comparison  laminate 
may  be  computed  from  the  general  results  of  Walpole  (1969),  specialized  to  the  case  of  isotropic, 
incompressible  phases.  We  obtain 


(4.3) 


1 
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where  the  overbars  on  the  moduli  denote  volume  averages  ie.g.,  =  and  where 

5=1 

Tj,  andr^  are  the  three  transversely  isotropic  invariants  of  the  applied  stress  tensor  a 

corresponding  to  the  three  independent  modes  for  an  incompressible,  transversely  isotropic  material 
(see  Appendix  I  and  Figure  1).  They  are  the  transverse  shear  stress,  the  deviatoric  stress,  and  the 
longitudinal  shear  stress,  respectively.  We  note  that  the  three  independent  modes  for  a  general 
incompressible,  transversely  isotropic  material  reduce  to  two  independent  modes  for  an 

incompressible  laminated  composite  (since  the  transverse  and  deviatoric  modes  have  the  same 
effective  response).  Note  further  that,  because  of  the  identity  from  the  section  on 

incompressible  materials  in  Appendix  I,  we  are  able  to  rewrite  the  first  term  in  brackets  in  (4.3)  in  the 
form 

With  expression  (4.3)  for  we  can  now  return  to  the  computation  of  U,  implied  by  (4.1). 
In  this  connection,  we  find  that  the  following  identity,  proved  in  Appendix  II,  is  useful  in  reducing 
the  number  of  optimizations,  namely. 


(4.4) 


1  .  I 

— =  nun  i 

n  (u<')  ,aisO  I 


where  the  (constant)  optimization  variables  (f  =  1, ...,  n)  are  required  to  satisfy  the  constraint 
5)  =  0.  Then,  substituting  (4.3),  together  with  (4.4),  into  (4.1)  leads  to  the  result 


(4.5) 

where 

(4.6) 


+ff. 


We  note  that,  by  definition,  the  functions 


are  concave  in  and 


similarly  the  variables  are  convex  in  q)^'\  Therefore,  by  the  Saddle  Point  Theorem  (Rockafellar 
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1970;  Corollary  37.3.1),  we  are  allowed  to  interchange  the  order  of  the  maximum  and  the  minimum 
in  (4.5).  Further,  it  follows  from  (2.8)  (assuming  convexity  of  where 
V,  =  1^)  that 

(4.7)  W" 

Therefore,  we  conclude  from  (4.5)  that 


(4.8) 


c/(S)  = 


min 

,  S>s0 


where  the  variables  are  given  by  relations  (4.6).  Evidently,  this  form  is  much  simpler  than  the 
original  form  given  by  (4.1)  and  (4.2);  it  involves  an  n-dimensional  constrained  optimization  in  place 
of  the  2n-dimensional  optimization  problem  implied  by  the  original  form.  However,  the  linear 
constraint  SJ  =  0  for  the  /i  optimization  variables  (r  =  1,...,  n)  can  be  embedded  into  the 
optimization  problem  (4.8)  by  letting  the  n-th  variable  be  expressed  in  terms  of  the  other  n  -  1 
variables  {s  =  1,...,  n  -  1)  via 


(4.9) 

With  this  modification,  the  problem  (4.8)  reduces  to  an  (n  -  l)-dimensional  optimization  problem 
over  the  unconstrained  variables  {s  -  1,...,  n  -  1).  For  instance,  for  the  case  of  a  two-phase 
laminated  composite,  the  problem  (4.8)  reduces  to  the  one-dimensional  optimization  problem 


(4.10)  (7(a)  = 


nun 


which  is  expressed  in  terms  of  one  (unconstrained)  optimization  variable  to.  Here,  we  have  made  the 
following  identifications,  =  c^^’o)  and 

Finally,  we  remark  that  simple  expressions  for  the  effective  stress/strain  relations  of  the 
nonlinear  transversely  isotropic  laminated  composite  may  be  obtained  by  means  of  the  results  of 
Appendix  HI.  These  relations  may  be  written  in  terms  of  the  incompressible,  transversely  isotropic 
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invariants  of  the  average  strain  tensor  e,  namely,  the  transverse  shear  strain  y^,  the  deviatoric  shear 
strain  and  the  longitudinal  shear  strain  7„  •  These  strain  invariants  are  defined  in  Appendix  I,  and 
are  completely  analogous  to  the  corresponding  (incompressible)  transversely  isotropic  invariants  of 
the  average  stress.  Thus,  with  the  help  of  relations  (in.7),  we  may  write 


(4.11) 


1  dy/^'^ 


L  --=1 


\2  1 


L  r=l 


y«  = 


^(D  \  ) 


L  /-'i 


s 

2 


where  and  where  the  variables  are  the  optimized  values  of  the  from  (4.8). 

We  note  that  for  the  nonlinear  laminated  composite,  there  is  full  coupling  between  all  the  distortional 
(shear)  modes.  This  is  different  from  the  situation  for  the  corresponding  linear  laminated  composite 
[see  (4.3)],  where  all  three  modes  are  uncoupled.  As  we  will  see  in  the  ensuing  discussions,  this 
inter-mode  coupling  is  one  of  the  intrinsic  features  of  laminated  (and  other  anisotropic)  nonlinear 
composites. 


5.  The  compressible  laminated  composite 

With  the  insight  gained  in  the  previous  section,  we  attempt  in  the  present  section  to  obtain 
corresponding  results  for  n-phase  laminated  composites  with  nonlinear,  isotropic,  compressible 
phases.  In  this  case,  we  can  apply  the  results  (3.3)  and  (3.4)  of  §  3  directly;  we  only  require  an 
expression  for  the  effective  energy  function  of  the  linear-elastic  laminate  with  isotropic, 

compressible  phases  in  prescribed  volume  fractions.  This  energy  function  may  be  computed  directly 
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from  the  results  of  Walpole  (1969)  for  the  transversely  isotropic  moduli  of  linear-elastic  laminated 
composites.  The  final  result  may  be  written  in  the  form 


(5.1)  Uoia)  =  U,(a)  +  U,ia), 

where 


with  T}^  =  9Kji^f[3K^+4^^),  and  where  and  t„  are  the  four  transversely  isotropic 

invariants  (up  to  quadratic  in  order)  of  the  applied  stress  a.  They  denote,  respectively,  the  in-plane 
hydrostatic  stress,  the  normal  tensile  stress,  the  transverse  shear  stress  and  the  longitudinal  shear 
stress  (see  Appendix  I  and  Figure  1).  The  reason  behind  the  above  splitting  of  U„  lies  in  the 

similarity  between  the  first  part  of  (5.1),  as  given  by  (5.2),  for  the  distortional  (shear)  modes  of  the 
compressible  laminate  and  relation  (4.3)  for  the  incompressible  composite  (with  replaced  by 

T^).  Thus,  it  follows  immediately  that 


where  the  to,  are  the  corresponding  optimization  variables,  and  they  are  subject  to  the  constraint* 
5),  =  0 .  The  second  part  is  more  complicated,  but  it  can  be  shown  by  straightforward  computation 
that,  if  may  be  represented  in  the  form 

(5.5)  (7,(5)  =  ™„  I -  (1  -  + 1(1  -  |. 

where  the  optimization  variables  are  also  subject  of  the  constraint  =  0. 
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By  putting  together  relations  (5.4)  and  (5.5),  we  arrive  at  the  following  expression  for  the 


linear  comparison  laminate 


(5.6) 


<»!:*.  s-=o 


1=1  r=l 


where 

(5.7) 


We  note  that  this  result  is  reminiscent  of  the  type  of  result  tliat  one  would  expect  to  arise  directly 
from  the  principle  of  minimum  complementary  energy.  That  this  result  is  indeed  directly  obtainable 
from  the  principle  of  minimum  complementary  energy  is  demonstrated  in  Appendix  IV. 

Then,  following  a  procedure  similar  to  the  one  followed  in  the  development  of  expression 

(4,8)  for  the  effective  energy  function  of  the  nonlinear,  incompressible  laminated  composite,  but 
making  use  of  (3.3)  and  (3.4),  we  arrive  at  the  following  expression  for  the  effective  energy  function 
of  the  nonlinear  laminated  composite 


(5.8)  fy(o)=  min  ■ 

5>.=5«=0 

where  and  are  given  by  (5.7).  Here,  we  have  made  use  of  the  Saddle  Point  theorem 
allowing  the  interchange  in  the  order  of  the  minimum  over  the  0^'^  variables  with  the 
maximum  over  the  comparison  moduli  and  We  note  that  form  (5.8)  fo:^  the  effective 

energy  function  U  of  the  nonlinear  laminated  composite  is  a  direct  generalization  of  lorm  (5.6)  for 
the  effective  energy  function  0^  of  the  linear  comparison  laminate.  There  are,  of  course,  sound 

physical  reasons  for  this  similarity,  and  these  arc  discussed  in  Appendix  FV,  where  it  is  shown  that 
the  pii'cedure  used  in  the  derivation  of  (5.6)  for  the  linear  laminate  can  be  extended  to  give  an 

alternative  derivation  of  (5.8)  for  the  nonlinear  laminate.  In  Appendix  FV,  we  also  show  that  an 
alternative  form  of  (5.8)  is  possible,  which  is  not  subject  to  the  ^  0  restriction,  although  we 
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should  emphasize  that  the  above  form  is  valid  in  the  limit  as  ->  0  (it  is  just  not  valid  in  a 
pointwise  sense  at  =  0,  because  the  optimisng  variables  become  unbounded  in  that  limit). 

We  note  that  while  the  distortional  and  dilatational  modes  are  not  coupled  in  the  linear 
laminated  composite  [although  the  dilational  modes  are  coupled  among  themselves,  as  (5.3)  shows], 
the  same  is  not  true  of  the  nonlinear  laminated  composite,  as  relations  (5.8)  and  (5.7)i  clearly 
demonstrate. 

The  new  representation  for  the  effective  energy  function  of  a  nonlinear  laminated  composite 
U  can  be  seen  to  involve  only  a  2n-dimensional  optimization  problem  with  two  linear  constraints. 
This  is  major  reduction  in  order  compared  with  the  original  expressions  (3.3)  and  (3.4)  involving  a 
4«-dimensional  optimization  problem.  However,  as  noted  in  the  previous  section,  further  reductions 
are  possible  [to  a  2(n  -  l)-dimensional  optimization  problem]  by  embedding  the  constraints  directly 
into  the  optimization  problem  (5.8).  For  example,  for  the  case  of  a  two  phase  composite,  we  obtain  a 
result  involving  only  a  two-dimensional  optimization  problem  prescribed  in  terms  of  the  variables 


(U,,  Q)„  via 
(5.9) 


U{k5)=  minji 

a»,,<u.  i 


where  and  are  given  by  relations  (5.7)  with  tuj'’ = 

and 

Finally,  we  remark  that  simple  expressions  for  the  effective  stress/strain  relations  of  the 
transversely  isotropic  laminated  composite  may  be  obtained  by  means  of  the  results  of  Appendix  III. 

These  may  be  written  in  terms  of  the  transversely  isotropic  invariants  of  the  average  strain  tensor  e, 
the  in-plane  hydrostatic  strain  e^,  the  normal  tensile  strain  £„,  the  transverse  shear  strain  and  the 

longitudinal  shear  strain  f, .  These  are  defined  in  Appendix  I,  and  are  completely  analogous  to  the 
corresponding  transversely  isotropic  invariants  of  the  average  stress.  Thus,  with  the  help  of  relations 
(ni.8),  we  may  write 
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/■-I 


y>.= 


3x1'^ 


where  f*'*  =  and  where  a)l'\  are  the  optimized  values  of  (o^'\ 

cu^'*  from  (5.8).  We  note  that  in  this  form,  the  coupling  between  the  distortional  and  dilatational 
modes  in  the  nonlinear  material  become  evident  since  depend  on  a,,  and,  conversely, 

e„,  depend  on  f^. 


6.  Application  to  laminated  composites  with  power-law  constitutive  behavior 

In  this  section,  we  specialize  the  results  of  §  4  and  §  5  for  three  classes  of  laminated 
composites.  The  first  subsection  deals  with  the  case  of  an  incompressible  laminated  material  made  up 
of  layers  of  a  phase  with  "linear  plus  power-hardening"  constitutive  behavior,  reinforced  with  stiffer 
layers  of  a  linear-elastic  material.  In  the  study  of  these  incompressible  laminates,  we  will  emphasize 
the  coupling  between  different  distortional  loading  modes  arising  as  a  consequence  of  nonlinearity 
and  anisotropy  in  the  laminates.  The  second  subsection  is  dedicated  to  the  study  of  a  compressible, 
aluminum/alumina  laminate,  and  the  understanding  of  the  dilatational  modes  is  emphasized  in  this 
case.  The  third  subsection  deals  with  an  incompressible  laminated  composite  made  up  of  two 
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rigid/perfectlv  plastic  phases  with  different  yield  stresses;  it  is  interesting  to  note  that,  in  this  special 
case,  completely  explicit  results  are  obtained  for  the  effective  yield  function  of  the  laminate. 


6.1.  Incompressible  laminated  composites 

In  this  subsection,  we  consider  an  incompressible,  two-phase  laminated  composite 
characterized  by  the  following  constitutive  laws  for  the  two  isotropic  phases.  Phase  #1  is  governed 
by  "linear  plus  power-hardening"  constitutive  behavior  described  by  the  energy-density  function 

(6.1) 

Jo 

where 


(6.2) 


F<"(s)  =  e, 


^  a  ^ 

y 


Here  H  is  the  unit  step  function  (equal  to  0  when  i  <  (J^,  and  to  1  otherwise),  and  e„,  are  strain, 

stress  normalization  factors  such  that  cr^/£„  with  denoting  the  shear  modulus  of  phase 

#1.  Then,  the  function  represents  the  uniaxial  stress/strain  relation  of  phase  #1  under  simple 

tension  loading  conditions.  Thus,  the  behavior  of  phase  #1  is  linear  when  the  uniaxial  stress  is  lower 
than  the  yield  stress,  a^,  and  is  linear  plus  power-hardening  for  stresses  larger  than  cr^.  The  factor 

V3  in  (6.1)  is  needed  in  order  to  fit  the  isotropic  stress  invariant  T,  to  the  uniaxial  case.  Phase  #2  is 

linear  and  governed  by  the  quadratic  energy-density  function 


(6.3) 


where  is  the  shear  modulus  of  the  phase. 

With  the  above  constitutive  behavior  for  the  two  phases  (#  1  and  2),  which  are  prescribed  in 
volume  fractions  (1  -  c‘^^)  and  c^^\  respectively,  the  effective  energy-density  function  of  the 
incompressible  laminated  composite  may  be  expressed  in  dimensionless  form  via  the  relation 


(6.4) 


Ll  Ll  Zl 

TV  \  T  T  T  /r  1/'*^ 

*^0*0 


r  T  T  (T  /X 

o  o  o  o  * 


19 


where  the  specific  form  of  the  function  G  is  deternuned  fi-om  (4.10),  and  so 

that  T„/y„  Then,  the  relations  between  the  three  (incompressible)  transversely  isotropic 

stress  invariants  and  the  corresponding  strain  invariants  may  be  computed  from  (4.11).  These 
relations  are  presented  in  Figures  2  through  5  for  the  following  values  of  the  four  parameters 
appearing  in  (6.4): 

(2) 

■^  =  1,  ^  =  5,  n  =  3andc<'^=0.2. 

We  recall  that  there  are  only  two  independent  modes  for  the  incompressible  laminated  composites; 
they  are  the  longitudinal  shear  stress  T,  and  the  following  combination  of  the  other  two  shear  modes 

(j.c.  the  transverse  and  deviatoric  shear  modes,  respectively).  For  simplicity,  we  will  refer 

to  this  combination  of  the  two  modes  as  t  and  to  the  corresponding  combination  of  the  strain 
modes,  +  yj ,  as  y.  Thus,  it  suffices  to  consider  the  relations  among  the  stress  modes  t,,  t 
and  the  strain  modes  y, ,  y  in  order  to  have  a  complete  description  of  the  constitutive  behavior  of  the 
incompressible  laminate.  In  order  to  highlight  the  effect  of  nonlinearity,  results  are  included  in 
Figures  2  to  5  in  the  form  of  short-dashed  curves  for  a  linear  laminated  composite  with  the  same 
shear  moduli  as  the  nonlinear  laminate.  Thus,  the  phases  of  this  linear  reference  laminate  are  similar 
to  those  of  the  nonlinear  one  with  the  only  difference  that  in  phase  #1  (7^  =  » . 

Figtire  2  shows  a  plot  of  the  longitudinal  shear  stress  versus  the  longitudinal  shear  strain 
y„  for  three  different  values  of  t  (t  /  =  0, 2, 5).  We  observe  that  when  there  is  no  preloading  of 

the  laminate  ( T  /  =  0),  the  behavior  of  the  stress/strain  curve  of  the  nonlinear  laminate  is  initially 

the  same  as  that  of  the  reference  linear  laminate  (short-dash  line)  until  phase  #  1  yields.  After 
yielding,  the  two  curves  diverge  with  the  nonlinear  phase  controlling  the  behavior  for  large 
longitudinal  shear  stresses.  That  this  should  be  so  is  seen  from  the  fact  that  shear  parallel  to  the 
layers  should  be  controlled  by  the  less  stiff  phase  (in  this  case,  the  nonlinear  phase).  The  effect  of 
increasing  T  is  to  saturate  the  linear  range  of  phase  #  1,  making  the  effective  stress/strain  curve  of 
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the  laminate  be  controlled  by  the  nonlinear  phase  even  for  small  values  of  the  longitudinal  shear 
stresses  t„. 

Figure  3  shows  a  plot  of  ?  versus  the  longitudinal  shear  strain  for  two  different  values  of 
the  longitudinal  shear  stress  =  0.5, 1),  and  serves  to  emphasize  the  coupling  between  the  two 
shear  modes.  Thus,  a  small  preload  in  the  form  of  a  longitudinal  shear  stress  applied  to  the 
nonlinear  laminate  can  lead  to  large  increases  in  the  longitudinal  strain  as  the  other  shear  stress 
mode  f  is  increased;  in  fact,  the  growth  is  unbounded  and  can  be  shown  to  be  proportional  to 

Figure  4  shows  a  plot  of  the  shear  stress  T  versus  the  shear  strain  y  for  three  different 
values  of  longitudinal  shear  stress  =0,0.5, 1).  We  observe  that  when  there  is  no 

preloading  of  the  laminate  =0),  the  behavior  of  the  effective  stress/strain  curve  of  the 

nonlinear  laminate  is  initially  the  same  as  that  of  the  reference  linear  laminate  (short-dash  line)  until 
phase  #  1  yields.  After  yielding,  however,  the  two  curves  diverge  with  the  linear  phase  controlling 
the  behavior  for  large  shear  stresses  t.  In  fact,  it  can  be  demonstrated  that  the  slope  of  the 
stress/strain  curve  in  question  reaches  an  asymptotic  value  of  (corresponding  to  a  linear 

Voigt  estimate  with  0)  as  the  shear  stress  f  becomes  large.  Evidently,  the  weaker  nonlinear 

phase  is  acting  as  if  it  was  not  present  for  large  enough  t.  The  effect  of  increasing  is  then  to 

saturate  the  linear  range  of  phase  #  1 ,  reducing  the  effect  of  the  nonlinear  phase  on  the  effective 
stress/strain  curve  of  the  composite  (the  laminate  behaves  almost  linearly  with  modulus  for 

sufficiently  large  preload  ?„).  That  the  nonlinear  laminate  should  be  controlled  by  the  stiffer  linear 
phase  for  large  magnitudes  of  the  transverse  shear  stress  (and  fixed  longitudinal  shear  stress  t,) 
is  easy  to  visualize,  but  that  exactly  the  same  behavior  should  be  observed  for  the  deviatoric  mode  f ^ 
(the  other  component  of  ?)  is  perhaps  less  intuitive.  The  reason,  however,  is  related  to  the  Poisson 
effect.  Thus,  for  example,  if  the  laminate  is  compressed  along  the  normal  direction  (which  may  seem 
to  be  controlled  by  the  less  stiff  nonlinear  phase),  tensile  strains  are  set  up  in  the  plane  of  the  layers. 
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which  must  be  continuous  across  the  phases,  thus  providing  the  required  stiffening  effect  in  the 
normal  direction  (because  the  linear  phase  controls  the  in-plane  behavior  of  the  laminate). 

Figure  5  shows  the  relation  between  the  longitudinal  shear  stress  mode  and  the  strain 
mode  y  for  different  values  of  the  shear  stress  t  (t  /  =  0.5, 2).  We  observe  that  whUe  there  is 

significant  coupling  between  the  two  modes  (by  comparison  with  the  linear  reference  laminate),  the 
coupling  is  not  as  significant  as  in  Figure  3.  Thus,  the  shear  strain  f  reaches  a  maximum  level  for  a 
given  shear  preload  T  as  the  longitudinal  shear  is  increased.  This  is  because  the  nonlinear  phase 
is  dominated  by  the  linear  phase  in  this  mode  of  deformation  as  observed  previously  in  connection 
with  Figure  4.  The  effect  of  increasing  preload  t  is  to  increase  (in  both  absolute  and  relative  terms) 
the  increments  in  the  shear  strain  y  with  increasing  shear  stress  f„. 

6.2.  The  aluminum/ alumina  laminated  composite 

In  this  subsection,  we  demonstrate  the  behavior  of  a  nonlinear,  compressible,  laminated 
composite  made  up  of  aluminum  layers  reinforced  with  layers  of  alumina.  Aluminum  is  a  ductile 
material  with  uniaxial  stress/strain  curves  that  can  be  approximated  by  a  "linear-plus-power"  law 
with  hardening  exponent  n  varying  between  4.2  and  5.8.  Thus,  we  will  assume  the  following  form 
for  the  energy-density  function  of  the  aluminum  layere  (phase  #  1) 

(6.5)  F^%)  ds  +  ol , 

where  is  the  same  as  in  (6.2),  and  thus  the  only  difference  between  (6.1)  and  (6.5)  is  the 
compressibility  of  aluminum  accounted  for  in  (6.5)  through  the  bulk  modulus  Alumina  (phase  # 
2)  is  a  brittle  material  that  behaves  in  a  linear  fashion  up  to  the  point  of  failure.  Its  energy-density 
function  is  represented  by 

(6-6) 

where  and  denote  the  shear  and  bulk  moduli  of  the  alumina,  respectively. 
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With  this  choice  of  and  (for  the  behaviors  of  the  two  phases),  the  effective 
energy-density  function  of  the  composite  can  be  represented  in  dimensionless  form  via 


(6.7) 


m  =  v"',  \. 


T^oYo 


<j  a  T  X  u 


(1) 


where  G  is  obtained  from  (5.9),  x^  =  Yo  -  '^oIYo  =  and  are 

the  (dimensionless)  Poisson's  ratios  of  the  two  phases  defined  by 


6K^'U2fi^'^ 


In  the  results  to  follow,  we  will  make  the  following  choices  (which  are  representative  of  the 
aluminum/alumina  composite)  for  the  material  parameters  in  (6.7): 


v<‘^=0.35,  v<^>=0.25  andAi  =  5. 


The  results  are  presented  in  Figures  6  through  9  in  terms  of  plots  of  the  four  transversely  isotropic 
stress  modes  versus  the  corresponding  strain  modes  for  three  different  values  of  the  volume  fraction 
of  alumina  (0.1,  0.25  and  0.5). 

Figure  6  shows  a  plot  of  the  in-plane  hydrostatic  stress  versus  the  corresponding 
hydrostatic  strain  e^,,  when  all  other  stress  modes  vanish,  for  the  three  values  of  the  volume  fraction 
of  alumina  .  We  observe  that  the  laminate  has  a  linear  range  with  effective  modulus  Iff  (recall 
that  T]  =  9Kfj./{3K+4iii))  up  to  yielding  of  the  aluminum  phase.  However,  the  laminate  behaves 
almost  linearly  even  after  yielding  with  modulus  approaching  for  large  values  of  This 


behavior  is  expected  on  physical  grounds  due  to  the  fact  that  the  stiffer  material  (alumina)  should 
dominate  the  behavior  in  tension  (compression)  parallel  to  the  layers.  The  effect  of  increasing  volume 
fractions  of  alumina  is  of  course  to  stiffen  the  effective  behavior  of  the  composite. 

Figure  7  shows  a  plot  of  the  normal  tensile  stress  versus  the  corresponding  tensile  strain 
e„,  when  all  other  stress  modes  vanish,  for  the  three  values  of  the  volume  fraction  of  alumina 
The  structure  of  the  plots  is  very  similar  to  that  of  Figure  6;  however,  the  effective  moduli  are 
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different.  Before  phase  #  1  reaches  yielding,  the  laminate  has  uniaxial  modulus  given  by  the 
expression 


while  after  yielding  of  phase  #  1,  the  modulus  for  large  stress  o;  is  reduced  to  the  level 

In  this  case,  it  is  not  evident  that  the  linear  phase  should  govern  the  effective  behavior  of  the  laminate 
for  large  stresses.  The  reason,  however,  is  the  same  as  discussed  in  the  previous  subsection  in 
connection  with  Figure  4:  continuity  of  the  tangential  strains  across  the  interfaces  together  with  the 
Poisson  effect. 

Figure  8  shows  the  corresponding  plots  for  the  transverse  shear  stress  versus  the 
transverse  shear  strain  when  no  other  stress  modes  are  present,  for  the  three  values  of  the 
volume  fraction  of  alumina.  In  this  case,  the  results  are  similar  to  those  of  Figure  6  for  clear  physical 
reasons:  the  stiffer  phase  controls  the  behavior  of  the  laminate  under  transverse  shear  loading. 

Figure  9  shows  plots  of  the  longitudinal  shear  stress  T„  versus  the  corresponding  strain 
mode  ,  with  no  other  stress  modes  present.  The  behavior  in  this  case  is  dramatically  different,  as 
the  study  of  the  corresponding  case  for  the  incompressible  laminate  demonstrated  earlier  (Figure  2). 
Thus,  after  an  initial  linear  range  before  yielding  of  phase  #  1,  the  weaker  nonlinear  phase  governs 
the  effective  behavior  of  the  laminate.  In  contrast  to  the  other  three  modes,  we  observe  that  the 
dependence  on  the  volume  fraction  of  alumina  is  fairly  weak,  so  that  the  three  curves  (corresponding 
to  different  values  of  are  quite  close  to  each  other. 

Qearly,  a  study  of  the  inter-relations  between  the  different  modes  would  be  required  to  have 
a  complete  picture  of  the  effective  behavior  of  the  nonlinear  compressible  laminate.  However,  the 
behavior  of  these  inter-modal  relations  is  similar  to  those  already  explored  for  the  incompressible 
laminate.  Thus,  the  inter-modes  relations  that  involve  the  longitudinal  shear  strain  y,  are  of  the  form 
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of  the  relations  presented  in  Figure  3  while  all  other  inter-modal  stress/strain  relations  are  in  the  form 
of  Figure  5. 


6.3.  The  rigidiperfectly  plastic  laminated  composite 

In  this  subsection,  we  consider  the  case  of  an  incompressible  laminated  composite  made  up 
of  two  rigid/perfectly  plastic  phases  with  yield  stresses  and  chosen  such  that  \  in 

given  volume  fractions  and  c^^\  The  behavior  of  the  phases  may  then  be  characterized  in  terms 
of  the  convex  energy-density  functions 


(6.8) 


Vr(')(T,)  = 


fo,  T,  <  ri^K 


(  r  =1,  2),  where  T,  denotes  the  effective  shear  stress.  These  energy  functions  may  be  obtained 
directly  from  pure  power-law  energy  functions  of  the  form 


(6.9) 
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y  \  *)  n+l  "  ‘ 
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in  the  limit  as  Further,  these  energy  functions  defrne  "yield  functions"  for  the  phase 

materials  that  may  be  described  in  the  usual  way  via 
(6.10)  0^^i<a)=T,-T<;^=O. 


Here,  we  will  proceed  formally  and  make  use  of  expression  (4.10)  to  determine  an 
expression  for  the  effective  energy  function  of  the  laminated  material  U ,  from  which  we  will  be  able 
to  determine  a  yield  function  for  the  laminated  composite  d>.  For  a  rigorous  treatment  of 
homogeneization  theory  for  rigid/perfectly  plastic  composites,  and  in  particular  for  a  discussion 
concerning  the  validity  of  the  normality  condition  for  the  effective  yield  function  of  the  composite, 

the  reader  is  referred  to  Suquet  (1985).  Because  of  incompressibility  and  transverse  isotropy,  we 
will  find  that  such  a  yield  function  may  be  represented  as  a  curve  in  the  (t„,t) -space,  where 

^  +  Xj  .  Thus,  application  of  (4.10)  leads  to  the  expression 

(6.11)  f7(a)  =  ntin{c(Tp,T„a))}, 
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where 


(6.12) 

and  where  = 


0,  ^  and  ^ 

oo,  or 

+  T^  and  =  ^(l  +  +  f^  . 


The  above  optimization  problem  for  O)  then  reduces  to  determining  all  possible 
combinations  of  T  and  T„  for  which  U  =  0,  which  in  turn  defines  the  yield  function  for  the 

composite  <I>.  First,  we  note  that,  independent  of  T  and  O),  U  can  only  vanish  if 

(6.13)  T„-zl‘><0, 

for  otherwise  ^  Thus,  inequality  (6.13)  is  a  necessary  condition  for  U  to  vanish. 

However,  the  condition  (6.13)  is  not  sufficient  to  ensure  that  U  vanishes  since  the  condition 
^(2)  <  ^2)  jjg  violated.  Thus,  assuming  that  condition  (6.13)  is  satisfied,  we  ask  the  question  of 
whether  there  are  values  of  to,  depending  on  f  and  T„,  such  that  conditions  and 

^  are  satisfied  simultaneously.  The  answer  is  affirmative,  provided  that  t  and  (for  given 
volume  fi^ctions  c^‘*  and  c^^*)  satisfy  the  condition 

(6.14) 


Thus,  conditions  (6. 13)  and  (6.14)  define  an  effective  yield  function  for  the  composite,  d>  =  0,  such 
that 


(6.15) 


We  note  that  when  rf  ’  =  \  the  expression  above  reduces  to  the  von  Mises  yield  criterion. 

Plots  of  the  yield  surfaces  in  the  (T„,T)-space  of  applied  stresses  are  given  in  Figures  10  and 
11.  Figure  10  shows  the  exact  yield  surface  <I>  for  the  choice  of  parameters,  =2  and 

=0.5.  The  isotropic  Reuss  and  Voigt  (also  known  as  Bishop-Hill  estimate)  bounds  for  the  yield 
surfaces  are  also  given  for  comparison.  We  note  that  the  exact  yield  surface  <I>  is  close  to  the  Voigt 
upper  bounding  surface  for  low  values  of  the  longitudinal  shear  stress  (?„  <  ^  t[*^),  and  close  to 
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the  Reuss  lower  bounding  surface  <!>,  for  low  values  of  t.  We  also  include  in  Figure  10  an  estimate 
for  the  yield  surface,  4>;,,  which  is  based  on  the  approximation  of  Hill  (1951)  for  slightly  anisotropic 
materials.  This  approximate  yield  surface  is  given  by 


(6.16) 


r 


+  1 
T+  ,  ~1» 


f  Wf 


and  we  note  that  it  amounts  to  an  elliptic  interpolation  between  the  Voigt  and  Reuss  yield  functions. 

We  observe  that  for  the  largely  anisotropic  case  depicted  in  Figure  10  (  rf  V  =  2  and 

=0.5),  Hill's  elliptic  approximation  severely  underestimates  the  ultimate  yield  strenght  of  the 
laminated  composite  for  combined  longitudinal  and  transverse  loading.  Figure  11  shows  plots  of  the 
exact  yield  surfaces  (continuous  lines)  and  Hill’s  approximate  yield  surfaces  (dashed  lines)  for  a 
laminated  composite  with  slight  anisotropy  (tf  V =1-25)  and  three  values  of  (0.1,  0.5  and 
0.9).  For  all  values  of  the  exact  yield  criterion  bounds  a  larger  region  of  the  (T„,T)-plane  than 


the  Hill  approximate  criterion,  and  the  two  curves  are  only  in  good  agreement  for  small  volume 
fractions  of  the  stronger  phase  (c^^^  =  0.1). 


6.  Closure 

In  this  paper,  we  have  described  the  application  of  a  new  variational  method,  developed  by 
Ponte  Castaneda  (1991a,  c),  to  a  laminated  composite  material  with  elastoplastic  phases  in  prescribed 
volume  fractions.  It  constitutes  one  of  the  first  applications  of  the  method  to  composite  materials  with 
anisotropic  symmetries  (see  also  Ponte  Castaneda  1991c  for  corresponding  results  for  fiber- 
reinforced  materials).  Because  of  the  simplicity  of  the  laminated  microstructure,  allowing  for  the 
determination  of  the  exact  effective  properties  of  the  laminated  composite,  this  work  is  of  interest  not 
only  on  account  of  the  practical  significance  of  laminated  composites  themselves,  but  also  because  it 
provides  a  case  study  where  the  power  of  the  method  can  be  evaluated.  Additionally,  the  results  of 
this  paper  suggest  that  when  dealing  with  strongly  anisotropic  materials,  it  is  not  sufficient  to 
consider  the  behavior  of  the  composite  under  special  loading  conditions,  since  the  behavior  of  the 
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composite  may  be  dramatically  different  under  a  different  type  of  loading  conditions.  Thus,  we 
found  that  nonlinearity  highlighted  the  differences  between  the  constitutive  response  of  the  laminated 
composite  under  transverse  and  longitudinal  shear  loading.  Further,  the  study  of  the  nonlinear 
laminated  composites  also  underlined  the  signiHcant  coupling  that  may  arise  between  different 
loading  modes  (which  does  not  take  place  for  linear  composites  with  similar  microstructures).  Thus, 
it  was  found  that  a  small  fixed  preload  of  the  laminate  in  the  longitudinal  direction  leads  to  continued 
increase  of  the  longitudinal  shear  strain  as  the  transverse  shear  stress  (for  example)  is  incremented.  It 
is  anticipated  that  the  features  uncovered  by  the  present  analysis  of  nonlinear  laminated  composites 
will  also  be  important  in  other  types  of  nonlinear  composites  with  anisotropic  symmetries,  such  as 
the  practically  important  class  of  fiber-reinforced  materials. 
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Appendix  I:  On  the  characterization  of  transversely  isotropic  materials 


The  purpose  of  this  appendix  is  to  gather  some  results  relevant  to  the  analysis  of  linear-elastic 
materials  with  transversely  isotropic  symmetry.  These  results  are  used  extensively  throughout  the 
body  of  the  paper  in  the  development  of  effective  stress/strain  relations  for  nonlinear  laminated 
composites,  which  are  a  special  class  of  transversely  isotropic  materials.  The  connection  between  the 
linear  and  nonlinear  laminated  composites  is  provided  by  the  new  variational  principle  in  terms  of  the 
exact  relation  (3.3).  The  emphasis  of  this  section  is  on  representations  for  the  transversely  isotropic 
invariants  of  the  stress  and  strain  tensor.  The  reason  is  that  nonlinear  transversely  isotropic  materials 
are  most  efficiently  characterized  in  terms  of  energy-density  functions  depending  on  these  invariants. 

I.l  Isotropic  invariants 

As  is  well-known,  there  are  three  isotropic  invariants  for  a  symmetric,  second-order  tensor. 
However,  only  two  of  these — those  that  are  of  quadratic  order,  or  less — ^are  relevant  to  linear-elastic 
behavior.  These  invariants  may  be  expressed  (see,  for  example,  Walpole  1981)  in  terms  of  two 
fourth-order  projection  tensors  J  and  K,  such  that  I  =  J  +  K,  JJ  =  J,  KK  =  K  and  JK  =  0.  Their 
Cartesian  components  are  given  by 
(1.1)  =  + 

where  5,^  is  the  Kronecker  delta  symbol.  Then,  in  terms  of  these  projection  tensors,  we  define  two 

isotropic  invariants  of  the  stress  tensor  via 

(1.2)  ~  1 

called  the  hydrostatic  (mean)  stress,  and  the  effective  shear  stress,  respectively.  We  also  define  the 
hydrostatic  strain  e„,  and  the  effective  shear  strain  y,  by  relations  completely  analogous  to  G.2). 

It  is  important  to  note  that  the  elasticity  tensor  L  of  an  Isotropic,  linear-elastic  material  admits 
a  spectral  decomposition 

(1.3)  L  =  3jfJ  +  2//K, 
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where  J  and  K  play  the  role  of  the  eigenprojections,  and  the  bulk  and  shear  moduli  of  the  material,  k 
and  fi,  are  the  corresponding  eigenvalues.  As  we  will  see  next,  the  situation  for  transversely 
isotropic  materials  is  different 

II.2  Transversely  isotropic  invariants 

There  are  in  general  five  transversely  isotropic  invariants  of  a  symmetric,  second-ordCT  tensor 
(Spencer  1971).  However,  only  four  of  these  invariants  are  linear,  or  quadratic,  in  order.  They  may 
be  represented  in  terms  of  the  four  projections  tensors  (see  Walpole  1981)  E*^',  and  E*^', 

satisfying  the  relations  E^^’e''*  =  E^'^;  E^'^E'**  =  0,  p  q;  and  E^'J  +  E^'J  +  E^'^  +  E''’^  =  I.  The 
components  of  these  four  projections  tensors  are  given  respectively  by 
(I't) 

where  =  n^n.  and  =  5^  -  n^nj,  with  n  denoting  the  axis  of  transverse  isotropy.  Then,  the  four 
transversely  isotropic  invariants  of  the  stress  tensor  o  may  be  expressed  in  the  forms 


{7K+<y22)}. 

{  ^33  }  ’ 

(1.5) 

1 

+ 

{(afj-KT^)}, 

which  correspond  physically  to  the  in-plane  hydrostatic  stress,  the  normal  tensile  stress,  the  (in¬ 
plane)  transverse  shear  stress,  and  the  (anti-plane)  longitudinal  shear  stress  (given  in  brackets  are  the 

corresponding  representations  for  a  choice  of  n  aligned  with  the  3-direction).  Analogous  relations 
apply  for  the  transversely  isotropic  invariants  of  the  strain  tensor  e,  denoted  respectively  f,, 

and  y,.  We  also  note  for  latter  reference  that  the  following  two  relations  hold  between  the 

transversely  isotropic  invariants  of  (1.5)  and  the  isotropic  invariants  of  (1.2),  namely. 


33 


(1.6) 


<T,  =  i(2<T,  +  <T.).  =  + 

Contrary  to  the  situation  for  isotropic  materials,  the  above  four  projection  tensors  are  not  the 
eigentensors  of  the  spectral  decomposition  of  an  arbitrary  transversely  isotropic  material  (Mehrabadi 
and  Cowin  1990).  Such  eigentensors  would  unfortunately  involve  the  material  moduli.  Therefore,  it 
is  necessary  to  introduce  (see  Walpole  1981)  two  other  tensors,  that  are  not  projections,  and 
with  components 

(1.7) 

Then,  the  elasticity  tensor  L  of  an  arbitrary  transversely  isotropic  material  may  be  expressed  in  terms 
of  these  six  tensors.  It  is  worth  mentioning  that  the  above  tensors  satisfy  the  relation 

(1.8)  J  =  |E'‘>  +  jE'"’  +  i(E>*l  +  E>®’), 

and  that  we  can  additionally  define  for  latter  reference  an  additional  tensor  E'  such  that 

(1.9)  E'  =  E'^>  +  E^^^-K, 

This  last  tensor  is  a  projection  tensor,  which  is  orthogonal  to  E'^’  and  E'^l 

Finally,  we  remark  that  the  energy  density  function  of  a  transversely  isotropic,  linear-elastic 
materials  may  be  represented  in  the  form 

a  10) 

Then,  the  relation  between  the  transversely  isotropic  stress  and  strain  invariants  is  given  by 


(I.ll) 


1  dy/ 
Ida/ 


and 


r»  = 


1  By 

2dx/ 


1.3  Incompressible,  transversely  isotropic  invariants 

For  incompressible,  transversely  isotropic  materials,  it  suffices  to  consider  the  three 
invariants  of  order  less  than  quadratic  on  the  space  of  traceless,  symmetric,  second-order  tensors. 
These  may  be  obtained  in  terms  of  the  three  orthogonal  projection  tensors  E*^',  E'^'  and  E',  defined 
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in  the  previous  subsection.  Thus,  the  incompressible,  transversely  isotropic  invariants  of  the  stress 
tensor  a  are  T^,  t„,  and  the  deviatoric  shear  stress 


(1.12) 


). 


corresponding  to  the  three  above  projections,  respectively.  We  note  further  that  from  (1.6)2  have 

the  following  identity  relating  the  effective  shear  stress  and  the  incompressible,  transversely  isotropic 
invariants,  +  t^.  The  corresponding  strain  invariants  are,  of  course,  denoted  y^,  y„ 

and  /rf. 


Finally,  we  note  that  the  elasticity  tensor  L  of  an  incompressible,  transversely  isotropic, 
linear-elastic  material  admits  a  spectral  decomposition  of  the  form 
(1. 13)  L  =  +  2n,E', 

where  fij,  are  the  three  shear  moduli  that  suffice  to  characterize  the  behavior  of  such  a 

material  (Lipton,  1991a). 


Appendix  II:  A  useful  identity 


In  this  appendix,  we  demonstrate  the  following  identity 


(n.i) 


1= 


which  is  used  repeatedly  in  the  body  of  the  paper. 
We  begin  by  letting  g  be  the  function  defined  by 


(11.2) 


r=l  ^ 


where  the  variables  (r  =  1,...,  n)  are  subject  to  the  constraint  co  =  1,  and  where  the  variables 
a*"’  >0  (r  =  1,...,  n). 
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The  choice  of  the  set,  =  a^'^a ,  satisfies  the  constraint  and  is  such  that  =  l/a . 

Consider  next  a  second,  arbitrary  set,  distinct  from  die  first  set,  (r  =  1,...,  n),  such  that  ^  =  1, 
and  let  6^'^^  =  -  (O^'^ .  Then,  substitution  of  this  second  set  into  (II.2)  leads  to 


(n.2) 


where  we  have  used  the  fact  that  0=0.  Hence,  identity  (I.l)  is  demonstrated.  In  the  body  of  the 
paper,  we  replace  by  (1  -  (o^'\  with  an  appropriate  modification  for  the  constraint 


Appendix  HI:  A  simplified  expression  for  the  effective  stress/strain  relations 


Consider  the  following  form  for  the  effective  energy  function  (4.8)  of  the  incompressible 
laminated  composite 


where  -  (0^‘^'f  +  T* .  As  shown  in  the  body  of  the  paper,  we  can  eliminate  the 

constraint  ©  =  0  by  letting 


(in.2) 


C  r=l 


and  rewriting  (111.1)  in  terms  of  the  /i  - 1  optimizations  variables  ©^'^  (r  =  l,...,n  - 1)  via 


ail.3)  t/(o)=  nun  + 

where  the  variables  (s  =  l,...,n-l)  are  the  same  as  before,  but  on  the  other  hand 


^(0©(0  ^  ^2  Then,  the  n-l  optimization  conditions  of  (III.3)  are 


given  by  the  relations 


au.4) 


=  0,  (r  =  !,...,« -1). 
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If  we  now  denote  the  optimal  variables  satisfying  (111.4),  by  (r  =  l,..,n  - 1),  the 
effective  energy  function  of  the  incompressible  laminate  may  then  be  written  in  the  form 

(in.5)  U{o)  = 

/a  ' 

where  j  and  = 

It  follows  that  effective  stress/strain  relations  of  the  laminated  composite  may  be  computed 
from  the  relations 


(III.6)  (f”)  (>+ Jol-''®'’)  +••• 

We  note  that  each  of  the  terms  in  the  last  summation  of  (III.6)  is  identical  to  zero  by  virtue  of  the 
optimizations  conditions  (in.4).  Thus,  in  the  computation  of  the  effective  stress/strain  relations,  we 
may  regard  the  optimizations  variables  as  constants  as  far  as  derivatives  with  respect  to  a  are 
concerned,  to  obtain  the  final  result 


(III.7) 


5 (f ''((1  -  % . 


where  is  defined  via  the  relation  (in.2)  in  terms  of  the  other  (r  =  1,...,  n  -  1). 

It  can  be  shown  that  an  analogous  result  may  be  obtained  for  the  nonlinear  compressible 
composite  with  effective  energy  function  U  given  by  (5.8).  In  fact,  we  may  write  the  effective 
stress/strain  relations  for  the  nonlinear  compressible  laminate  in  the  form 


(1II.8) 
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where  and  0^1^^  are  evaluated  from  (5.7)  at  the  optimal  values  of  (ol'’\  and  denoted  by 
and  respectively.  Here,  the  derivatives  with  respect  to  the  average  stress  o  are  evaluated 
with  d)l'\  and  fixed. 

Appendix  IV:  More  on  the  compressible  laminated  composites 

Having  obtained  expressions  (5.6)  and  (5.8)  for  the  effective  energy  functions  of  the  linear 
and  nonlinear  laminated  composites,  and  we  note  that  the  form  of  these  expressions  is 

reminiscent  of  the  type  of  result  that  one  would  expect  from  direct  utilization  of  the  principle  of 
minimum  complementary  energy  (2.4).  In  this  appendix,  we  briefly  show  that  results  (5.6)  and  (5.8) 
for  and  U ,  respectively,  can  indeed  be  alternatively  obtained  directly  from  the  principle  of 

minimum  complementary  energy.  It  is  important  to  emphasize,  however,  that  while  the  derivations 
given  in  the  body  of  the  paper  result  from  straightforward  mechanical  computations,  the  present 
derivations  based  on  the  principle  of  minimum  complementary  energy  rely  on  the  physics  of  the 
problem,  and  were  motivated  by  the  prior  derivations  we  would  not  have  discovered  them  had 
they  not  been  derived  by  the  prior  analysis).  Also,  the  case  of  a  laminated  composite  is  a  very  special 
case  of  a  composite  material;  in  general,  we  do  not  expect  that  other  types  of  results,  such  as  bounds 
for  the  effective  properties  of  nonlinear  isotropic  composites,  obtained  by  means  of  the  new 
variational  principle  (e.g.  Ponte  Castaneda  1991a,  b,  c),  may  also  be  obtained  directly  from  the 
minimum  complementary-energy  principle. 

We  begin  with  the  derivation  of  the  linear  result  (5.6).  We  have  already  mentioned  that  the 
stress  field  within  the  laminated  composite  is  piecewise  constant,  i.e.,  of  the  form 

o  =  •  n)a^'’^ ,  where  corresponds  to  the  constant  stress  field  in  phase  r.  The  problem 

r«l 

then  reduces  to  that  of  finding  these  unknown  phase  stresses  together  with  the  corresponding 
constant  strain  fields  (related  to  the  stresses  by  the  phase  constitutive  relations),  and  satisfying 
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the  conditions  of  continuity  of  the  traction  stresses  and  tangential  strains  across  the  interfaces 
between  the  phases,  as  well  as  the  averaging  conditions  stated  in  §  2. 

In  this  connection,  the  interior  and  exterior  projection  operators  of  Hill  (1972,  1983) 
F  =  and  E  =  E^^'  +  E^^^  (refer  to  Appendix  I),  respectively,  turn  out  to  be  useful  because 

they  allow  the  decomposition  of  any  symmetric,  second-order  tensor  into  its  tangential  (interior)  and 
traction  (exterior)  components  (with  reference  to  a  boundary  with  normal  n).  Thus,  the  tangential 
components  of  the  strain  (which  must  be  continuous  across  inteiphase  boundaries  on  the  laminated 
composite)  are  given  by  Fe,  and,  correspondingly,  the  traction  components  of  the  stress  (which 
must  also  be  continuous  across  the  interphase  boundaries)  are  given  by  Ea.  Alternatively,  we  may 
state  that  E^^^o,  E^^'a  and  E^‘^e,  E^^'e  must  also  be  continuous  across  such  boundaries. 

Next,  we  apply  the  above  results  to  the  laminated  composite,  for  which  the  interfacial 
boundaries  are  all  perpendicular  to  a  fixed  vector  n.  Since  the  traction  stresses  must  be  continuous 
from  phase  to  phase,  we  have  that 

(IV.  1)  EfV'’  =  Ef%,  and  EfV'^  =  Ef"^o, 

where  we  have  additionally  made  use  of  the  average  stress  condition  given  in  §  2.  We  continue  by 
noting  that  for  an  isotropic  material  (as  are  all  the  phases  in  the  our  laminate), 

within  each  linear  phase,  and  therefore  for  an  isotropic  phase  the  E^^'  projection  of  the  stress  tensor 
must  have  the  same  direction  in  all  phases.  Thus,  applying  the  averaging  condition  for  the  stresses, 
we  arrive  at 

(IV.2)  =  (I  -  )E^'>a , 

where  the  variables  must  satisfy  the  condition  that  57,  =0.  Additionally,  since  the  E^‘^- 
projection  is  one-dimensional,  it  follows  that  the  E*‘^ -projections  of  the  stress  tensor  must  also  be 
parallel  from  phase  to  phase.  Therefore,  applying  the  averaging  condition  for  the  stresses,  we  have 
that 

(IV.3)  =  (1  - 
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where  the  variables  (O^'^  must  satisfy  the  condition  that  (0^=0.  We  note,  however,  that  if  =  0 
(or,  equivalently,  if  o'p  =  0),  the  above  result  does  not  hold,  because  in  this  case  the  corresponding 

projections  of  the  stress  in  the  phases  need  not  vanish  (only  their  average  needs  to  vanish). 

Applying  the  results  of  Appendix  I  [in  particular,  (1.6)],  we  conclude  that  the  isotropic 
invariants  of  the  stress  tensor  within  each  phase  and  (on  which  the  energy-density  functions 

of  each  isotropic  phase  depend)  are  precisely  those  given  by  relations  (5.7).  Therefore,  it  follows 
from  the  principle  of  minimum  complementary  energy — ^by  minimizing  over  the  set  of  admissible 
stresses  (i.e.,  over  the  optimizing  variables  and  subject  to  the  constraints  =  0  and 
=  0) — that  the  effective  energy  function  of  the  linear  composite  is  indeed  given  by  expression 

(5.6). 

For  the  nonlinear  laminated  composite,  we  observe  that  the  same  analysis  given  above  would 
also  work,  leading  to  expression  (5.8)  for  U.  The  only  modification  that  is  required  in  this  analysis 
is  that  for  a  nonlinear  isotropic  phase  (say  phase  r),  the  relation  would  not 

hold,  but  it  can  be  easily  shown  that  for  the  nonlinear  isotropic  material  of  the  type  considered  in  this 
work,  the  conclusion  (rV.2)  would  still  hold,  and  hence  the  final  form  for  U  would  be  the  same  as 
that  for  the  linear  laminated  composite  . 

We  conclude  this  appendix  by  stating  an  alternative  form  of  (5.6)  and  (5.8)  that  works  even 
when  =  0).  This  is  accomplished  by  redefining  the  optimizing  variables  in  terms  of  the  new 

variables 


— _ E_^2! _ 1 


(IV.4) 

where  now  we  need  to  have  that  -  a,  0.  In  terms  of  the  new  variables  (the  variables 
do  not  change),  relation  (5.8)  is  expressed  in  the  form 


av.5) 


U{o)= 
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where 


(IV.6)  f w  =  ^(l  -  rl  +  +  i[(a,  -  CT,)(l  -  and  ai*^  =  <t,  +  |(<t,  -  <t,)(i  -  • 


Note  that  when  =  cr,  =  0,  we  are  guaranteed  that  =  cf['^  =0  in  each  phase,  and  then  both 
forms  are  equally  valid. 
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List  of  figure  captions 
Fig.  1.  A  two-phase  laminated  material. 

Fig.  2.  The  relations  between  the  longitudinal  shear  stress  and  strain  of  the  incompressible, 

nonlinear  laminate  (continuous  lines),  and  the  reference  linear  laminate  (short-dashed  lines)  for  three 
different  values  of  the  other  stress  mode  ?;  f/T^  =  0,  t/Tq  =  2,  and  t/Tq  =  5. 

Fig.  3.  The  inter-relations  between  the  shear  stress  T  and  the  longitudinal  shear  strain  of  the 

incompressible,  nonlinear  laminate  (continuous  lines),  and  the  reference  linear  laminate 
(short-dashed  lines)  for  two  different  values  of  the  longitudinal  shear  stress  t„;  t„/To  =0.5,  and 
VTo  =  1. 


Fig  4.  The  reladons  between  the  shear  stress  t  and  the  conresponding  shear  strain  y  of  the 
incompressible,  nonlinear  laminate  (continuous  lines),  and  the  reference  linear  laminate 
(short-dashed  lines)  for  three  different  values  of  the  longitudinal  shear  stress  t„/To  =0, 
T„/To=0.5,and  t„/To=2. 

Fig.  5.  The  inter-relations  between  the  longitudinal  shear  stress  f„  and  the  shear  strain  f  of  the 

incompressible,  nonlinear  laminate  (continuous  lines),  and  the  reference  linear  laminate 
(short-dashed  lines)  for  two  different  values  of  the  stress  mode  t;  t/Tq  =  0.5,  and  t/Tq  =  1 . 

Fig.  6.  The  relations  between  the  in-plane  hydrostatic  stress  and  strain  of  the  compressible 
laminate  (where  no  other  stress  modes  are  present)  for  three  different  values  of  the  volume  fraction 
of  the  linear  phase, 

Fig.  7.  The  relations  between  the  normal  tensile  stress  and  strain  of  the  compressible  laminate 
(where  no  other  stress  modes  are  present)  for  three  different  values  of  the  volume  fraction  of  the 
linear  phase, 

Fig.  8.  The  relations  between  the  transverse  shear  stress  and  strain  of  the  compressible 
laminate  (where  no  other  stress  modes  are  present)  for  three  different  values  of  the  volume  fraction 
of  the  linear  phase,  . 
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Fig.  9.  The  relations  between  the  longitudinal  shear  stress  and  strain  7,  of  the  compressible 
laminate  (where  no  other  stress  modes  are  present)  for  three  different  values  of  the  volume  fraction 
of  the  linear  phase,  . 

Fig.  10.  The  exact  yield  surface  of  the  laminated  composite,  the  anisotropic  elliptic  estimate  of  Hill 
and  the  Voigt  and  Reuss  isotropic  estimates  plotted  on  the  (t„,t) -plane  for  =2  and 

c^^^=0.5. 

Fig.  11.  Tlte  exact  yield  surface  of  the  laminated  composite  (continuous  lines)  and  the  quadratic 
approximation  of  Hill  (dashed  Unes)  plotted  on  the  (t„t) -plane  for  rf  V rf  ^  =  1.25  and  three  values 
of  (0.1,  0.5  and  0.9). 
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